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CHAPTER I 
INTRODUCTION 


The text which follows constitutes the final report on work performed at the Sperry Rand 
Research Center (SRRC) under contract NAS 9-9705 with the NASA Manned Spacecraft Center, 
Houston, Texas, The objectives of this program were twofold. First, it was to consider opti- 
mum organization and component technology for a strapdown transformation computer based upon 
the assumption of gyroscope and accelerometer inputs denoting components of angular and veloc- 
ity change measured about an orthogonal axis triad. Second, it was to treat optimum methods 
of implementing a reliable attitude reference or control system when the gyroscope and accel- 
erometer inputs were redundant and configured about the six normals to faces of a dodecahedron 

The work which is reported here divides into four major activites, each corresponding to 
a chapter. These will be described briefly. The first part of the activity, described in 
Chapter 2, is a study of means whereby a strapdown computer such as the one constructed for 
laboratory use under a previous contract may best be organized and implemented in reliable, 
flight-ready form. The second activity, covered in Chapter 3, is an investigation of optimum 
means for performing attitude reference computations when redundant dodecahedron sensor arrays 
are utilized. The third part of the study, described in Chapter 4, is a study of an attitude 
control system employing redundant control moment gyros as actuators. This part of the work 
was performed at the Sperry Flight Systems Division, under subcontract to SRRC. The final 
part of the study, Chapter 5, contains the results of further computer simulations of the trun 
cation-error-free direction-cosine computation used in SIMU, methods and results are described 
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CHAPTER 2 

CONFIGURATION AND COMPONENT TECHNOLOGY FOR A STRAPDOWN COMPUTER 

In a previous contract a strapdown coordinate-conversion computer was designed and con- 
structed which used a truncation-error-free algorithm to compute direction cosines relating a 
rotating-axis system to one which was mertially stabilized. The organization of this com- 
puter was based upon a building block concept in which most blocks were digital accumulators 
with overflow detection and ternary input gating. In the present study, where the goal is to 
determine optimum system organization and component technology for such a computer, several 

t 

structures, including the accumulator block, were studied. 

Before describing the system, the underlying reasons for the particular choice of imple- 
mentation will be given. First of all, since discussion centers about a spaceborne computer, 
weight, volume and power must be minimized and reliability must be maximized. From the point 
of view of weight, volume and power, it is clear that minimality is provided by incorporating 
as high a level of circuit integration as is feasible within the constraints of present day 
technology. Studies have shown that maximizing the level of integration also provides the 
greatest reliability,'*' since circuit bonds and off-chip interconnections, which are the pri- 
mary sources of failure, are minimized. At the present time, reasonable device yields are 
obtained for integrated circuit chips with dimensions of 120 mils on each side. Interline 
spacings of 0.1 mil are readily obtainable. For metal-oxide-silicon (MOS) circuitry, active 
device areas are approximately 1.0 rail square. Allowing for metallization patterns and bond- 
ing pads, this means that an active device count of 2000 to 5000 is reasonable for MOS inte- 
grated circuits. Four -phase or ratioless MOS circuitry allows transfer rates from 2 to 5 MHz 
to be realized, and provides the highest circuit density presently available. For this reason, 
the preferred mechanization of the strapdown computer would use four -phase MOS LSI. 

As in most applications, cost is also of importance. Development and manufacturing costs 
for the strapdown computer can be reduced by limiting the number of different chip types 
required. For this reason and the reasons above, the approach which was followed in formu- 
lating the systems structures was one in which maximum use was made of the fewest possible 
different circuit types, while maintaining the highest level of integration consistent with 
present technology and while trying to minimize the number of interconnections between chips. 

The block diagram assumed for the triad system is shown m Fig. 2-1. It is essentially 
that of the SIMU computer, except that sections intended primarily for laboratory use have been 
eliminated. A set of three signals from orthogonal triad gyroscopic sensors provides the 
input to the direction cosine computation. These signals are assumed to be corrected for 
drift scale factor errors by equipment preceding the navigation computer. The gyro inputs are 
added in the buffer to triad rotational rate commands. These commands replace the earth-rate 
correction system in the SIMU computer and are assumed to be generated by the main guidance 
and navigation computer. They enter the strapdown computer referred to stable coordinates 
and are transformed to rotational coordinates by the "commanded rate transformation" unit. 

Triad strapdown accelerometer signals are transformed to stable coordinates by the "AV 
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FIG. 2-1 Block diagram of SIMU computer. 
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ti ansformation" unit, as in SIMU. System outputs arc the stabilized AV’s and the nine direc- 
tion cosines. 


Besides the configuration used in the SIMU computer, a new form has been developed for 
system implementation. This will be described first, and then methods for increasing system 
reliability will be discussed for both cases. 

The new method of organization is very similar in form to SIMU except that 16-bit digital 
resolvers are used for the direction cosine computation instead of the accumulators. A dia- 
gram of one of the three identical direction cosine units is shown in Fig. 2-2. Each direc- 
tion cosine resides in a triple of resolvers, the Y resolver storing the upper 16 bits, the 
R resolver storing the middle 16 bits, and the E resolver storing the lower 16 bits. 

The updating algorithm used is the original truncation-error-free method involving three 
passes. The operations for a 9^ pulse are specified by the matrix equation 


M 


X 



( 2 . 1 ) 


Similar operations hold for 9^ and 9^ . 

Each of the resolvers is equipped with its own serial adder operating in a closed loop 
with the storage register Thus, simultaneous accumulation by all resolvers is allowable. 

By judicious sequencing of operations, a three-to-one speedup in processing rate is obtained 
over the originally proposed SIMU computer. This sequencing is best clarified by means of the 
schedule shown in Fig. 2-3. Each term in the schedule signifies the updating of the contents 
stored in the specified resolver by adding to it the contents of another resolver, as deter- 
mined by the updating matrices. The superscripts indicate whether it is the X-, y-, or z-matrix 
equation, and the priming of the superscripts specifies the right, middle, or left matrix, 
respectively. Thus, E^ indicates the updating of resolver for the second pass of 9^ 
updating. This updating involves subtracting the contents of resolver Rg from the contents 
of resolver E 1 and storing the difference in E^ . It is clear that the updating for x, y 
and z pulses of any rank of resolvers (say, the E resolvers) requires just nine word-times. 
Since overlapping of cycles is allowable, a new set of gyro pulses may be entered every nine 
word-times. If the clock used is the same as in the SIMU computer, a processing rate speedup 
of two-to-one over SIMU is obtained The speedup over the originally proposed system is three- 
to-one. Besides the speedup in processing rate (or, conversely, the ability to use a lower 
clock rate), another advantage of this form of implementation is that it provides greater uni- 
formity of components throughout the system, since the resolvers used in the cosine computa- 
tion are identical to those used to provide the coordinate transformations elsewhere in the 
system. An appropriate cucuit embodiment would be the realization of one or more resolvers 
as a single integrated circuit. The approximate device count for a simplex realization of 
such a circuit is 3^10 transistors, so that such an embodiment is quite feasible using MOS/LSI 
technology. 
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FIG. 2-2 Resolver direction cosine unit. 



FIG. 2-3 Sequencing of all-resolver cosine computer. 
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Several techniques have been investigated for increasing reliability over the "simplex" 
SIMU computer. These will be described, although no revolutionary developments have emerged. 

The first technique makes use of the orthogonality properties of the direction cosine matrix 

T 

for self checking. If C is the cosine matrix and C is its transpose, then orthogonality 
implies that 

C • C T = C T • C = I (2.2) 

where I is the identity matrix. Carrying out the term-by-term multiplications, the twelve 
distinct equations in Table 2.1 result. 


TABLE 2.1 

Direction Cosine Orthogonality Conditions 


Eq. No. 

Equation 


2 2 2 

1 

C 11 + C 12 + C 13 - 1 = 0 

2 

C 11 C 21 ~ C 12 C 22 + C 13 C 23 = 0 

3 

C 11 C 31 + C 12 C 32 + C 13 C 33 = 0 


2 2 2 

4 

C 21 + C 22 T C 23 - 1 = 0 

3 

C 21 C 31 + C 22 C 32 + C 23 C 33 = 0 


2 2 2 

6 

C 31 + C 32 + C 33 - 1 = 0 


_2 J2. „2 

7 

C 11 + C 21 + C 31 - 1 = 0 

8 

C 11 C 12 + C 21 C 22 " C 31 C 32 = 0 

9 

C 11 C 13 + C 21 C 23 + C 31 C 33 = 0 


„2 „2 2 

10 

°12 + '"22 + c 32 “ 1 " U 

11 

C 12 C 13 + C 22 C 23 * C 32 C 33 = 0 


2 2 2 

12 

°13 + C 23 + C 33 ~ 1 ~ U 


The failure of some subset of the equations in Table 2.1 to be satisfied may be used as 
a form of parity check to determine faulty computation of the direction cosines. For instance, 
an error m C.^ alone is indicated by lack of satisfaction of Eqs, (2.1) and (2.7). Similar 
criteria may be derived from the covering relationships expressed by Table 2.2. At "X" at 
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the intersection of a particular column and row indicates that the cosine corresponding to that 
row is covered by the equation corresponding to the column. 

TABLE 2.2 

Direction Cosine Covering by Orthogonality 
Equations (see Fig. 2-1) 



If a single bit error is introduced into one of the direction cosines, say C^ , the error 
will propagate to cosines C^ 0 and C^g as a function of angular input increments received. 

No error, however, will propagate to the other cosines, since their computation "loops” do not 
intersect those of C^ . The dependence of C-^ . for instance, upon may be expressed 
by the partial differential equation pair below 



12 


fa. . 

ae ' ''ll 

Z 


(2.3) 


(2.4) 


Consider the inversion of a single bit m the representation of C- ^ and its effect upon 
C^2 • The bit inversion will be equivalent to adding some quantity S to C^ . Let the error 
in be and the error in C^ be ^ ’ Then 


*u 3E l 

90 + 39 -C 12 " E 2 


(2.5) 
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( 2 . 6 ) 


30 38 c n E 1 

z z 


with initial values 



Subtracting left and right equivalents from this set, the resulting error equations are 

(2.7) 

( 2 . 8 ) 



with the obvious solution 


E, - S cos 6 (2.9) 

1 z • 

E 0 = S sin 9 (2.10) 

2 z 

If the checking equations of Table 2.1 are solved at fixed intervals and involve the upper n 

bits of the cosines, then the maximum period between solutions such that a single bit error in 

one of the direction cosines can be detected before it affects the other cosines in its loop is 

a function of n and the maximum angular rate which the system can accommodate. If the worst 

case bit inversion (i.e., S=l) is assumed, the maximum period between solutions is set by the 

inequality below, where uj is the maximum angular velocity and T is the maximum period 

in3x max 

T £ — i- sm" 1 2‘ n . (2.11) 

max id 

max 

# 

After detecting the cosine which is in error, the problem of what to do about it still 
remains. If the angles were unchanging, it would be a simple matter to recompute the correct 
\nlues of the cosines using the relationships of Table 2.1. This, however, is not generally 
the case. One possibility is to retain in storage at some point the last complete set of cor- 
rect cosines computed before the error was detected. These values can be reinserted and the 
computations restarted. Although this method introduces an overall error if vehicle attitude 
has changed during the time since restart values were observed, the amount of error can be 
limited if the checking equations are run frequently enough. Another method which leads to no 
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error accumulation is to stop the direction cosine computations as soon as an error is detected 
and store the input angle increments in a buffer, maintaining the proper sequence among the 
three axes. The equations in Table 2.1 then can be used to restore the proper value to the 
faulty cosine, after which the stored input increments are utilized to update the cosine matrix 
to its proper value. This method will succeed only if the allowable updating rate of the 
cosine matrix is greater than the maximum allowable pulse input rate from the gyros. 

In the event that the error is caused by permanent failure of one of the components within 
the direction cosine computer, more than one-time correction will bo required. A possible 
solution to this situation is obtained by noting that the nine cosines are computed as three 
independent triads. If a spare triad computer is included and is equipped with the proper 
connections to allow its output to be switched to replace any of the three primary computer 
outputs, then upon detection of a failure in one of the three primary units the restart opera- 
tion can involve the spare unit, after which it is switched into the system in place of the 
faulty element. 

In the methods described above that make use of the orthogonality properties of the direc- 
tion cosine matrix for self correction, the nature of the computations and their rate is such 
that corrections can be handled by thegeneral purpose guidance computer associated with the 
system. It is possible, and also more straightforward, to include logical redundance within 
the direction cosine computer itself rather than using external means. In both the SIMU and 
all-resolver systems, numerous computational feedback loops are used. If it is desired that 
the addition of redundancy will effect correction of transient errors as well as permanent com- 
ponent failures, it is advisable to add the redundancy in such a way that error correction 
occurs within all feedback loops. One scheme for providing this sort of correction is shown 
in Fig. 2-4, which depicts a triply modular redundant resolver. Voting circuits are introduced 
in such a way that they intersect the tightest feedback paths and still provide correct outputs 
from all three of the resolver segments for the maximum number of component failures within the 
segment . 

It is difficult to supply all-inclusive arguments for the level at which redundancy should 
be introduced in the strapdown computer, but since the resolver is the element which seems best 
suited for use as the basic building block, larger or smaller elements not providing the same 
lead-count and circuit-type minimality, and since the packaging of three interconnected resol- 
vers with voting elements on a single substrate is presently within the state of the art, this 
configuration is recommended as the form to be used for implementation of the strapdown com- 
puter . 

REFERENCES 

1. "Failure Rate/Temperattirc Data for Umvac Military Computers," Univac Federal Systems Divi- 
sion, Publication PX-438B-2, July 1969. 
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CHAPTER 3 

CONVERSION METHODS FOR TRANSFORMATION FROM DODECAHEDRON-TO -TRIAD AXES 

3.1 DIRECT REALIZATION BY USE OF PSEUDO-INVERSE AND STATUS MATRICES 

3.1.1 Theoretical Analysis of Approach 

In this part of the report we describe a method of transforming the six-vector 


(3.1) 


of the dodecahedron reference system into the three-vector 

(3.2) 

of the orthogonal triad reference system. 

The geometric relationship between the two systems is given by 

m = Hb (3.3) 





- 11 - 



Since every 3 x 3 submatnx of H is nonsingular, any three out of the six equations of 
(3.3) can be used to determine the b components in terms of the corresponding three components 
of m . The nature of the problem, however, is such that the instruments providing the com- 
ponents of m are subject to errors and, as was pointed out by Gilmore, the solution of Eq. 
(3.3) for b which will best fit the geometric configuration of the well functioning instru- 
ments is given by 


b = (h T X 1I T A m (3.7) 

where X = diaq(X , X. ,X ,X.,X ,X.) is the "status" matrix, with -1 or 0 according to 
3 a b c d e f ' 1 J 

whether instrument i, i=a ,b,c,d,e,f , is error free or not. 

Of course, for this solution to be meaningful, at least three out of the six instruments 
has to be error free. There are 42 different combinations with at least three of the X^^ 
being equal to 1 . After having evaluated explicitly the matrix 


rl 1 = (H T X H j H T X 


(3.8) 


for all the proper 42 combinations of the X^ , it becomes evident that the total number of 
different nonzero entries of H* for the various combinations is only 25. In terms of S and 
C , these 25 entries are as shown in Table 3,1. 

TABLE 3.1 

Twenty Five Required Data-Word Matrix 


N 1 

= — ^2S 
10 

+ c ) 

N ll = 

* 

S + 4c) 

N 21 ~ 

2S 

; + c 

«2 

= To( 7S 

+ c) 

N 12 = 


3S 

< 

+ C) 

N 22 = 

K 

3S - C) 

N 3 

- i's + 

■ 5 V S 

3C) 

N 13 = 

3 

4 

S 


N 23 = 

* 

S + 3c) 

N 4 

tO 

\ 

II 

+ 7Cj 

N 14 = 

3 

4 

c 


N 24 = 


s - cj 


“Ik 5 - 

- 2Cj 

N 15 = 

1' 

f 

C S 

+ 3c) 

N 25 = 

K 

s.c) 


N 6 = §(3S - C) 
"t = To[ 2S ~ c ) 

% =i s 

»10 = - c ) 


16 


K 2s - c ) 

«„ * K S + 2c) 

N l8 = |)-S + 2C) 

= 1(2S + c) 


»20 - -S ♦ 2C 
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The 42 explicit evaluations of H 1 are given in Table 3.2 in terras of the indices of 
the above listed N.'s . For instance, the matrix for the combination anc * 

i\. c =l^=\ f =0 , which is given by 


N 19 

- N !9 

0 

0 

0 

0 


-N 22 

" N 22 

0 

0 

N 20 

0 

(3.9) 

_ N 18 

N 18 

0 

0 

0 

0 



In Table 3.2 ,0 stands for the number zero and the numbers shown represent the terras by 


listing only the subscript number. (±i). 

e.g. 

’ * N 12” 

: 12 . 

For the 

expression m Eq. 

(3.9), this 

yields 

r- 









19 

-19 

0 

0 

0 

0 



X ,X. ,X = 1 

-22 

-22 

0 

0 

20 

0 


(3.10) 

a b e 









18 

18 

0 

0 

0 

0 



3.1.2 Implementation of Pseudo-Inverse 

and 

Status 

Matrices 

Method 

Osina LSI 



Packaged Basic Elements 


The implementation of the transformation described m 3.1.1 essentially involves the use 
of a read-only memory and an arithmetic unit to provide the appropriate overflows to the 3-axis 
system’s input buffer logic. 


The following discussion will treat each section of the system in enough detail to make 
clear the functional interaction of the sections in achieving the required final results. 
There are four sections of the system: they are (1) the six input gyro-sync logic, (2) the 

read-only memory, (3) the arithmetic unit, and (43 the direction cosine computer's (SIMU) 
input buffer logic. 

The gyro-sync logic performs two operations. First, the six (ix, iy, ±z) asynchronous 
gyro inputs are synchronized to computer word time. Second, these six signals are encoded so 
as to be sent out as three octal variables (see Fig. 3-1). 

The read-only memory's address selection is a function of two sets of input variables. 
These are encoded gyro outputs described above and the set of variables ( 1 j, ,Xg and 
X^} that indicate which gyros are operating satisfactorily at any given time. These nine 
selection lines are used to select one of 252 memory locations. In each of these memory loca- 
tions is an address of a second read-only memory. The second read-only memory contains 25 
sixteen-bit words. 


To understand the need for this memory configuration refer to Table 3.2, which shows the 
manner in which the basic constants shown m Table 3.1 are used. The entries in Table 3.1 are 
a summary of the various products of a constant times ^ Sln a sad cos ce as discussed m 
Sec. 3.1.1. 
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FIG. 3-1 Incremental input encoding 
for the six-gyro system. 
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TABLE 3.2 
Status Matrices 




“ 19 

-19 

0 

0 

0 

0" 



'21 

0 

0 

0 

-23 

23“ 

W X c ~ 

1 

-23 

23 

21 

0 

0 

0 

X a' X e ,X f = 

i 

0 

0 

0 

0 

18 

18 



_18 

18 

0 

0 

0 

0_ 



_ 0 

0 

0 

0 

19 

-19. 



"19 

-19 

0 

0 

0 

0“ 



■ 0 

0 

18 

18 

0 

0“ 

X a’ X b’ X d = 

1 

23 

-23 

0 

-21 

0 

0 

X b ,X c ,X d = 

i 

0 

0 

19 

-19 

0 

0 



_18 

18 

0 

0 

0 

0. 



_ 0 

20 

22 

22 

0 

0_ 



" 19 

-19 

0 

0 

0 

0“ 



' 0 

18 

23 

0 

-19 

0“ 

\ ^ “ 
a’ fa* e 

1 

-22 

-22 

0 

0 

20 

0 

W X e = 

i 

0 

-19 

-18 

0 

23 

0 



_ 18 

18 

0 

0 

0 

0. 



_ 0 

23 

19 

0 

-18 

0_ 



~ 19 

-19 

0 

0 

0 

0" 



" 0 

-22 

19 

0 

0 

-18“ 

VV X f = 

1 

22 

22 

0 

0 

0 

20 

X b ,X c ,X f = 

i 

0 

18 

22 

0 

0 

19 



_ 18 

18 

0 

0 

0 

0- 



_ 0 

19 

18 

0 

0 

-22. 



“ 0 

0 

18 

18 

0 

0" 



- 0 

-22 

0 

19 

18 

0" 

a c d 

1 

0 

.0 

19 

-19 

0 

0 

X b ,X d’ X e = 

i 

0 

-18 

0 

-22 

19 

0 



.20 

0 

-22 

-22 

0 

0- 



_ 0 

19 

0 

18 

22 

0. 



"22 

0 

19 

0 

-18 

0“ 



' 0 

18 

0 

23 

0 

19" 

\ ,X = 

ace 

1 

-18 

0 

22 

0 

19 

0 

X b ,X d* X f " 

i 

0 

19 

0 

18 

0 

23 



_ 19 

0 

-18 

0 

22 

0- 



- 0 

23 

0 

19 

0 

18. 



--18 

0 

23 

0 

0 

-19" 



- 0 

-21 

0 

0 

23 

-23“ 

X a ,X c ,X f = 

1 

19 

0 

-18 

0 

0 

23 

X b’ X e’ X f = 

i 

0 

0 

0 

0 

18 

18 



_ 23 

0 

-19 

0 

0 

18_ 



- 0 

0 

0 

0 

19 

-19. 



"-18 

0 

0 

23 

19 

0" 



“ 0 

0 

18 

18 

0 

0" 

^ X , f A. — 

1 

-19 

0 

0 

18 

23 

0 

\ ,\ Jt X = 

i 

0 

0 

19 

-19 

0 

0 

a d’ e 








c* d* e 










- 23 

0 

0 

-19 

-18 

0_ 



- 0 

0 

-23 

23 

21 

0. 



~ 22 

0 

0 

19 

0 

18" 



■ 0 

0 

18 

18 

0 

0 " 

X a' X d’ X f ~ 

1 

18 

0 

0 

-22 

0 

19 

c d f 

i 

0 

0 

19 

-19 

0 

0 



L_ 19 

0 

0 

-18 

0 

-22. 



„ 0 

0 

23 

-23 

0 

-21. 


- 15 - 



TABLE 3.2 

Status Matrices (cent.) 
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TABLE 3.2 

Status Matrices (coin.) 
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In operation, the first read-only memory determines which matrix in Table 3 2 will be 
required to process the sensors which are operable (1--1). The second memory selects the 
words specified by this matrix. These words are stored as in Table 3.1. They are processed 
with the appropriate timing to allow each of the triad axes to be updated for any gyro outputs 
that have occurred. These words are used in 42 different sequences in performing the trans- 
formations in the arithmetic unit. There are three sets of memories as described above for 
processing from six to three axes, one for each of the triad axes. 


The arithmetic unit consists of an adder and five data registers used m parallel process- 
ing to transform to 6 9x, 69y and 69 z (see Fig. 3-2). The sequence of events is as follows. 
Starting the updating to obtain a a9x input to the computer, the first specified word from 
the second read-only memory (ROM) (of the set of six words specified by the first ROM for the 
x-axis update)is loaded into the constant buffer. From the buffer the scaled increment is 
added to the previous remainder of A0x and the result is held in the general 68 buffer 
register. Once this is done the selection of the next constant for the x-axis update from the 
ROM can begin. The 69x register is loaded again to prepare for the next component of the 
A9x update. After the six cycles are completed the remainder will be in the A0x remainder 
register and an output will be sent to the gyro-buffer logic if the entire update produced 
either a positive or negative overflow. 


The gyro-buffer will recognize this overflow as the signal to produce a positive or nega- 
tive output pulse that will update the appropriate direction cosines. This sequence will be 
followed by the second set of ROM's and its arithmetic register to provide a A9y input to 
the direction cosine matrix, if the inputs are of sufficient magnitude to produce an overflov. 
The same sequence will take place for the final set of ROM’s to produce a A9z input pulse 
when required. 


17 - 



The above constitutes a complete update cycle and takes place within the system update 
time of 57.6 ps. The computer described above can function continuously at any input rate 
desired that does not exceed this pulse repetition rate, 

3.2 CONVERSION FROM DODECAHEDRON TO TRIAD BY SENSOR CORRECTION 

3.2.1 Description 

In this section a method for converting sensor signals from redundant dodecahedron to 
equivalent triad axes is described in which signals missing because of sensor failure are 
restored by synthesis front operational sensor outputs. The resulting six inertial signals 
are then used as the input to a fixed transformation whose output is the desired equivalent 
triad set. A block diagram of the system is shown in Fig. 3-3. It consists of a sensor 
failure detection unit, a correction unit, and a dodecahedron-to-tnad conversion unit. The 
detection unit is the parity equation integrator, which will be described in Sec. 3.2.2. The 
detection unit produces as its output a six-element binary vector, \ , each of whose elements, 
when equal to unity, affirms the operation of one of the sensors. The \ vector and the sen- 
sor outputs serve as inputs to the sensor correction unit. If a particular sensor is opera- 
tive, the output of the correction unit corresponding to that sensor is the sensor signal 
itself. If a sensor has failed, then the correction unit output corresponding to that sensor 
is an estimate of what the sensor signal should be, derived from the operative sensor signals 
It is clear that at least three sensors must be operative m order that such an estimate can 
be made. 

The six sensor signals (synthetic or actual) from the correction unit are used as inputs 
to the dodecahedron-to-tnad conversion unit. Since six sensor signals are always assumed to 
be present at the input of this unit, even when as many as three sensors have actually failed, 
a simple, fixed conversion algorithm may be used. The techniques employed and the system 
itself will now be described m greater detail. 

The method used for synthesizing failed-sensor signals is based upon solution of one of 
the parity equations for the missing signal in terms of the three valid signals. The parity 
equations, as defined by Gilmore, are reproduced here as Table '3.3. As an example, if sensor 
A has failed but sensors B, C , and D are still operating, then Eq. (1) in Table 3.3 may 
be used for estimating A , viz, 

A = B t Y(C + D) (3.11) 


where 


Y = StB) = ^Ir 1 = °- 618033 • (3.12) 

A is the estimate of the output that would be obtained from sensor A if it were operating 
and o is the angle between two intersecting normals to adjacent faces of a dodecahedron. A 
general equation may be written for sensor coirection : 
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1 = X X 1 + ^ 1 {^2 X 3 X 4^6 + X 2 X 3 X 4 X s) E A 
+ X 2 X 3 X 4 X 5 X 6 E B + X 2 X 3 X 5 X 6 E C 
+ X 2 X 3 X 4 X 5 E D + X 2 X 3 X 4 X 5 X 6 E E 

(3.13) 

+ ^2 X 3 X 4 X 5 X 6 E F + X 2 X 3 X 4 X 5 X 6 E G 
+ X 2 X 3 X 4 X 5 X 6 E 4 + ^3^ 4 ^5^ 6 E i 

+ X 2 X 4 X 5 X 6 E j} ; 

TABLE 3.3 

Failure Detection Parity Equations 


No. 

Instrument 

Equation 

1 

ABCD 

(A - 

B)c 

- 

(C 


D)s 

= 

0 

2 

ABCE 

(B + 

C)c 

- 

(A 


E)s 

= 

0 

3 

ABCF 

(C - 

A)c 


(B 

- 

F)s 

= 

0 

4 

ABDE 

CD - 

A)c 

- 1 - 

(B 


E)s 

= 

0 

5 

ABDF 

(B -r 

D)c 

- 

(A 

- 

F)s 

= 

0 

6 

ABCF 

(E - 

F)c 

- 

(A 

+ 

B)s 

= 

0 

7 

ACDE 

(D + 

E)c 

- 

(A 

- 

C)s 

= 

0 

8 

ACDF 

(F - 

C)o 

- 

(A 

+ 

D) s 

= 

0 

9 

ACEF 

(A + 

F)c 

- 

(C 

+ 

E)s 

= 

0 

10 

ADEF 

CE - 

A)c 

+ 

(D 

- 

F)s 

= 

0 

11 

BCDE 

(E - 

C)c 


(D 

- 

B)s 

= 

0 

12 

BCDF 

(F + 

D)c 


(B 

- 

C)s 

= 

0 

13 

BCEF 

(B - 

E)c 

+ 

(C 

4 . 

F)s 

= 

0 

14 

BDEF 

(B + 

F)c 

+ 

(D 

- 

E)s 

= 

0 

15 

CDEF 

(C - 

D)c 

- 

(E 


F)s 

= 

0 



(/ 5 4- 

5^ 








C = cos(ff) = 

V 10 

) 

~ 

O.U5065 













s = sin(oO = 

V 10 


** w 

0 . 

52574 




To interpret Eq. (3.13) let 1 represent the output of sensor 1 and let 1 represent 
the estimate of this quantity. The various X’s are the operational indicators of their 
respective sensors or, when barred, the complements of these quantities. Finally, the quanti- 
ties E A* e b , etc., represent the parity equations to be used in computing 1 . Thus, Eq. 
(3.13) serves to select which (if any) equation will be used to synthesize the sensor output, 
based upon the operational status of all sensors in the system. 
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Equation (3.13) can be made specific to any particular sensor by means of Table 3.4, 
which indicates the substitution of subscripts for correction of a given sensor. The synthesis 

of a failed-sensor output from one of the parity equations involves multiplying certain of the 

-1 „ 

operational outputs by either unity, Y or Y = 1+y and then summing these quantities. The 
selection of quantities to be summed to produce the desired synthetic sensor output is governed 
by the X vector. A resolver implementation of sensor output synthesis is shown in Fig 3-4 
for the general case. Six such resolvers are required, one for each dodecahedron axis. The 
resolver signals , either multiplied by Y or by unity, are gated into a resolver by 

sequencing signals T. and control signals C . The overflow output of this resolver is 
the dodecahedron sensor signal, either synthetic or real, to be used in the dodecahedron-to- 
tnad conversion. The control signals are derived from Eq. (3.13) as functions of the 

X vector and are given in general form by Table 3.5. Table 3.5 is made specific to a particu- 
lar sensor by means of Table 3.4. Implementation of the control functions is by means of a 
read-only memory whose inputs are the five X. and whose outputs are the C . This memory 
has a capacity of 16 words, each 15 bits long. Six identical memories are used, one for each 
of the dodecahedron sensors. Each memory is made specific to an axis by the combination and 
permutation of X's used as its input. These are specified by Table 3.4. 

TABLE 3.4 

Table for Substitution in General 1-Gyro Restoring Equation 


Sensor 

General Variable 

B 

X 2 

X 3 

X 4 

X 5 

X 6 

e a 

e b 

E c 

e d 

e e 

e f 

e g 

e h 

E I 

E J 

A 

X A 

X B 

x c 

X D 

X £ 

X F 

E 1 

E 2 

E 3 

E 4 

E 5 

E 6 

E 7 

E 8 

E 9 

E 10 

B 

X B 

X A 

X C 

X D 

X F 

X E 

E 1 

E 3 

E 2 

E 5 

E 4 

E 6 

E 12 

E 11 

E 13 

E 14 

C 

x c 

X D 

X E 

X F 

X A 

X 8 

E 15 

e t 

E 11 

E 0 

E 12 

E 1 

E 9 

E 13 

E 2 

E 3 

D 

X D 

x c 

X £ 

X F 

X B 

X A 

E 15 

E 11 

E 7 

E 12 

E 8 

E 1 

E 14 

E 10 

E 4 

E 5 

E 

X E 

X F 

X A 

X B 

x c 

X D 

E 6 

E 9 

E 10 

E 13 

E 14 

E 15 

E 2 

E 4 

E 7 

E 11 

F 

X F 

X E 

X A 

X B 

X D 

x c 

E 6 

e io 

E 9 

E 14 

E 13 

E 15 

C 5 

E 3 

E 8 

E 12 
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FIG. 3-4 Dodecahedron gyro correction resolver (6 required). 
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TABLE 3.5 

Table of Binary Constants for Dodecahedron Gyro Correction (see Fig. 3-4) 


Eq. 

X 2 

X 3 

*4 

X*. 

0 


C 21 

*2l 

C 22 

*31 

*31 

C 32 

C 41 

*41 

c 42 

C 51 

*51 

*52 

*61 

*61 

*62 

A 

1 

1 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 


1 

1 

1 

1 

0 

















1 

1 

1 

1 

1 
















B 

1 

1 

0 

1 

0 

1 

0 

1 

1 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

C 

1 

1 

0 

0 

1 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


1 

1 

1 

0 

1 
















D 

1 

0 

1 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 


1 

0 

1 

1 

1 
















E 

1 

0 

1 

0 

1 

1 

0 

1 

0 

0 

0 

1 

0 

1 

0 

0 

0 

1 

0 

0 

F 

1 

0 

0 

1 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

1 

1 

G 

0 

1 

1 

1 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

1 

0 

1 

0 

0 

0 

H 

0 

1 

1 

0 

1 

0 

0 

0 

1 

0 

1 

0 

1 

0 

0 

0 

0 

0 

1 

1 

I 

0 

1 

0 

1 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

1 

0 


1 

' 1 

0 

1 

1 
















J 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0, 

1 

1 

0 

0 

0 

0 

1 


0 

1 

1 

1 

1 

















When all six dodecahedron sensor signals are assumed present, the dodecahedron-to-triad 
conversion takes the form 
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(3.14) 


The implementation of the conversion is by means of three resolver circuits, each identi- 
cal to the one shown in Fig. 3-5. Whole-number representations of the constant half-sines and 

A 

half-cosines are gated into the resolver by synthetic sensor outputs , in accordance with 
Eq. (3.14) and by sequencing signals T. . The overflow output of each resolver is the 
respective triad signal. _ 23 - 





FIG. 3-5 Triad axis generator (3 required). 
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3.2.2 Parity Integration 


To implement the parity equation determination of gyro failure, a weighted integration 
technique will be adopted. The reasoning behind the use of integration is that quantization 
effects in both the gyros and the computer must be separated from significant gyro deviations. 
Weighting is included to mask out long-term drift due to roundoff and other effects within the 
parity computation itself. Consistent with other sections of the computer, an operational 
digital procedure will be used in the parity section of the system. 

The form of the parity conditions used closely resembles that shown by Gilmore. 1 That is, 
for the first partly equation 

= (s A - Sgjcos a - + Sg^sin cr (3.15) 

where is the output of gyro i and is a measure of ui . For perfect gyros, Q^O . 

Thus Q. represents an angular-rate eTror. Since pulsed gyros are assumed, where each' output 
pulse represents an increment of angle, the first parity equation may be written m angle form 
as 


dX l = ^l dt = ( d ®A " d ®B) C0S “ ( d6 c ^ d0 D > ln a (3.16) 

and similarly for the other equation. Temporal weighting and integration are introduced by 
defining the quantity P^ , 


P 

i 


J ax, 


0 


(3.17) 


or 



(3.18) 


If dX^dt-0 for t<0 then is seen to be the convolution of dX^/dt with e The 

differential form of the above integral equation may readily be shown to be 


dP 

dt 


-KP 

l 



or 


dP = -KP dt + dX 

ill 


(3.19) 


(3.20) 
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The final form implies a simple DDA solution, as shown in Fig. 3-6. 



FIG. 3-6 DDA solution. 

The implementation of this configuration using resolvers is shown ill Fig. 3-7. A complete 
solution to allow correction of two failures and detection of three simultaneously requires 15 
units of the type in Fig. 3-7, although the elements so designated may be shared among the 15 
systems. It is also possible to mechanize a single-failure correction scheme using just six 
of the units in Fig. 3-7 in which the dX j inputs are switched among the parity networks. The 
gyro failures (up to two) are indicated by the binary quantities X , i-A ,B,C,D,E,F . The 
value X =1 corresponds to the condition that gyro i is operative. The X. are binary 
functions of the parity threshold functions . Specifically, 

X = n (e ) (3.2D 

1 subset ** 
i 

The E subset for each X is indicated by Table 3.6. The E numbering corresponds to 
Jit J 

Table 2 in Gilmore, Detection of three or more failures is given by 

15 

X T = n E . (3.22) 

J=1 J 

The condition X^=0(X^,= 1) indicates that at least three gyros have failed, but the failures 
cannot be isolated. This is equivalent to complete system failure. 

The rationale behind the use of the weighted integral form of parity equation follows that 
discussed by Keene.“ The prevalent type of sensor failure is performance degradation, evidenced 
by either scale factor change or by excessive bias. Let g be the sensor output. Then 

g = (a + ®)<u + b (3.23) 

where a is the scale factor, e is the scale factor error, b is the bias, and ui is the 
measured quantity. Since integrating sensors are normally used, the output of interest is 

J gdt = (a + e) J uidt + bt . (3.24) 

Then for constant iu = iu . 
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FIG. 3-7 Parity integrator. 
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TABLE 3.6 
Subset for X. 



8 = I gdt = an t + ( e<o + b)t 
J a o o 


(3.25) 


implying a ramp error behavior. Since X. is a linear combination of the 8's , it too exhib- 
its ramp behavior under these conditions. Let X^Sgt . Then Eq. (3.19) becomes 


or 



(3.26) 



-Kt 


) 


(3.27) 


If M ns defined as the maximum allowable value of P , then threshold is exceeded for 


P > M 
i 
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or 


T (* - o'") » “ 


(3.28) 


Solving for the time at which the inequality is satisfied, 


t 


M 




(3.29) 


Thus, KM is established as the minimum error rate for parity failure, and the time until such 
failure is detected is adjusted by K and decreases with increasing Q o . 

3.2.3 Constant Rate Multipliers 

Although it has been suggested that standard resolvers be used to effect the multiplica- 
tion Df a sensor or other pulse-rate signal by a constant quantity, another interesting approach 
has been developed. The basis for this approach is given in Appendix A (Smooth Sequences) and 
stems from the realization that the pulse rate resulting from the multiplication of the source 
rate by a fixed quantity (less than unity) should preserve as well as possible the relative 
pulse density of the source. Although the use of resolver-type rate multipliers does in fact 
preserve relative densities optimally, several techniques have been developed which perform 
the same function with appreciable circuit savings. One of these well suited for applications 
in the dodecahedron sensor transformation computer will be described. Consider Fig. 3-8, which 
shows an interconnection of pulse counters and inhibit elements. The upper counter produces 
one output pulse for every a^ input pulses, and so forth. Each time a pulse counter produces 
an output pulse it inhibits one input pulse to the counter immediately above it. Analysis of 
this circuit shows that the pulse output frequency is equal to the pulse output frequency mul- 
tiplied by P/Q , where P and Q are relatively prime integers and 


P 

0 



(3.30) 


Equivalently, P output pulses are produced for every Q input pulses. The P output pulses 
so produced are synchronous with input pulses, but are spaced as uniformly over the input pulses 
as is possible. 

The representation of P/Q in Eq. (3.30) is as a simple continued fraction. The proper- 
ties of such fractions are well known. One property of particular interest here is the 
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representation of irrational numbers by their approximation m continued fraction form. 

Theory shows that all algebraic irrational numbers may be represented by an infinite continued 
fraction whose partial quotients (a.^ in Eq. (3.30)) are bounded. The theory shows further 
that the approximation of any rational or irrational number by a truncated simple continued 
fraction is optimal in the sense that the absolute error incurred is smallest over all rational 
approximations with denominator no larger than that of the continued-fraction-derived approxi- 
mation. 

A case of special interest m the dodecahedron-to-tnad conversion system is the multipli- 
cation of sensor output pulse rates by the quantities 


_ sin a _ /5 - 1 
Y cos a 2 

and 

V -1 = cos of _ /5 + 1 _ , + v 
Y sin or 2 y 


(3.31) 


(3.32) 


where 2® is the angle measured between any two dodecahedron normals. The expansion of Y as 
an infinite simple continued fraction is given in Eq. (3.33) 


Y = 


1 


+ 


1 

1 

1 + 

1 + 


_1 

1 

1 + . 


(3.33) 


This leads to a very simple pulse-rate multiplier of the type shown m Fig. 3-8, since all 
the a i are unity and are therefore just direct connections. The resolution of this multiplier 
is determined by the number of stages used. For a resolver-type multiplier the resolution is 
doubled with each stage that is added to the storage register, since this has the effect of add- 
ing one more bit to the number added. It is interesting to compare the resolution per stage of 
the resolver to that of the continued fraction multiplier. Clearly, the resolution of an n-bit 
resolver-type multiplier is l/2 n and is independent of the quantity by which the rate is mul- 
tiplied, The approximate resolution of the continued fraction multiplier can be calculated 
using known properties of continued fraction expansions. The n 1, convergent of a continued 
fraction is defined as the rational fraction obtained by truncating the continued fraction after 
the n** 1 partial quotient. If is the n t convergent of some continued fraction, then 

C, = P n /Q„ (3.34) 


and 



(3.35) 


C n 


_ (-D 


n -1 


'n -1 


Q Q 

n n- 


In Eq. (3.33) 


lim C = y 
n-°° n 


The odd-order convergents form a monotomc sequence that converges to Y from above and the 
even-order convergents form a monotomc sequence that converges to Y from below. The abso- 
lute error E = | C — Y I is therefore bounded by the difference between the n 1 "* 1 and n~l St 
n 1 n 

convergents. Letting this difference be , 


and 


D = C - C J = 7 -^ — 
n 1 n n- 1 1 Q n Q n _ 1 


1 D n-ll ^ C n-1 " C n-2^ Q ,C 


Q n-l°n -2 


(3.36) 


(3.37) 


\ 

but 


lim Q 
rr*“ 


n- 


1 


: Y Q„ 


Therefore, in the limit. 


D 

n 



(3.37) 


(3.38) 


Thus, each additional partial quotient included in the continued fraction approximation to Y 
increases the resolution by a factor of approximately Y = ( 3 +/5)/2 = 2 .618 , and as n 
increases the resolution approaches 1/Y~^ n = 1/(1 +Y) 2n . The resolution of resolver -type and 
continued-fraction rate multipliers is compared m Table 3.7. In the strapdown system a resolu- 
tion equivalent to a word length of 16 bits is generally required. It may be seen from Table 
3.7 that this resolution is obtained using only twelve stages of the continued fraction multi- 
plier. Another advantage of the continued fraction multiplier is that it is essentially a 
parallel system with no carry propagation delays. Therefore, slower electronic components may 
be used to attain the samespeeds as the equivalent resolver. 

A working model of a 12-stage continued fraction multiplier was constructed. The circuit 
diagram of this unit appears as Fig. 3-9. 
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All NAND gates TI 7400 
All NOR gates TI 7402 
All inverters TI 7404 

FIG. 3-9 Continued-fraction pulse rate multiplier 













TABLE 3.7 


Comparison of Resolution of Resolver 
and Continued Fraction Ratio Multipliers 


n 

Resolutior 

Cont. Fraction 
Resolution 

i 

2 

2 

2 

4 

6 

3 

8 

15 

4 

16 

40 

5 

32 

104 

6 

64 

273 

7 

128 

714 

8 

256 

1870 

9 

512 

4895 

10 

1024 

12816 

11 

2048 

33552 

12 

4096 

87841 

13 

8192 

229970 

14 

16384 

602070 . 

15 

32768 

1576239 

16 

65536 

4126648 

17 

131072 

10803704 
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CHAPTER 4 


INCREMENTAL COMPUTER CONFIGURATION 
FOR A REDUNDANT CMG CONTROL SYSTEM 

4.1 INTRODUCTION 

The objective of this study is to present an incremental computer structure 
for a space vehicle attitude control system that uses a set of six skewed, 
single-gimbal, control m oment gyros (CMG's) as actuators and a set of six rate 
gyros in a dodecahedron configuration as rate sensors. The computer is to em- 
ploy the redundancy in both the sensors and actuators to provide fail-operational 
performance. The CMG configuration of six skewed, single-gimbal gyros (exempli- 
fied by the Sperry 6-GAMS configuration) can be sized to provide three-axis 
control with failures in as many as three of its gyros, and the dodecahedron- 
configured rate gyro package can provide vehicle rate information from any three 
of its six rate gyros. 

If. 2 GENERAL DESCRIPTION OF THE CMG CONTROL SYSTEM 

Since the computations to be performed by the CMG control computer are di- 
rectly related to the specific configuration and steering law selected for the 
CMG’s and also to the rate sensor configuration, these items are described 
briefly in paragraphs 4.2.1, 4.2.2, and * 1 . 2 . 3 , respectively. Paragraph 4.2.4 
identifies the functions to be performed by the computer, including a discussion 
on three -variable versus six- variable input data to the steering law computer. 

A brief discussion on computer failure detection is presented in paragraph *f.2.J. 

If, 2.1 The CMG Conf iguration 

Previous studies conducted jointly by Sperry and Lockheed for the Air 
Force 1 have identified a family of skewed, single-gimbal CMG configurations, 
denoted as the GAMS family of CMG configurations, to he the most efficient in 
terms of weight, power, and volume in providing redundant three-axis attitude 
control. Specifically, the 6-GAMS configuration., illustrated in Figure 4.2-1, 
is particularly suitable where control is to be maintained with any two of the 
six CMG's not operating. Although the control computer in this study is directed 
principally toward this configuration, the results are expected to be adaptable 
to other configurations. 

The black arrows in the CMG configuration model shown in Figure 4-1 de- 
pict the angular momentum vectors j h^ j of the six gyros, all of approximately 
equal magnitude. These vectors are shown in their respective reference direc- 
tions in the figure. The notation used in describing the configuration angles 
is shown in Figure 4-2. All gimbal axes are tilted at an angle p from the 
x-axis, and the projections of the gimbal axes on the yz-plane are directed at 
angles y = j 0, 60, 120, 180, 240, 300 j degrees from the y-axis. The ginihal 
angles denoted by a = jc^, are measured from their respective reference 

directions in the counterclockwise direction about the gimbal axes as shown by 
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FIG. 4-1 Model of the Sperry 6-GAMS configuration. 


X 



FIG. 4-2 CMG configuration reference system. 
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the angular scales of each gyro in Figure 4-1. For a = jo, . of, the net 
angular momentum vector, H, of the CMG configuration is zero, assuming that all 
the individual momenta, jh^j, are of equal magnitude. In general, the projec- 
tions of H onto the vehicle principal axes are given hy equations (4.1), where 
S (• ) denotes sin ( ) and ^ denotes cos ( ). 


H_ = 2 
x i=l 


h i S p S a i 


H = 2 -h. fa C 
* i=l 1 V Yi a i 





(4.1) 


H = 
z 


2 h. (C C 
L=1 1 V Yi 



E 

Yi 



Equations (4.1) do not include the angular momentum due to the rotations of 
the gimbals and rotors abont the gimbal axes -which, for the purpose of steering 
law computations, may be assumed to be negligibly small in comparison with the 
angular momentum due to the rotation of the rotors about their spin axes. 


CMG systems provide attitude control by appropriate transfer of angular 
momentum between the CMG'S and the vehicle. The mathematical law hy which de- 
sired vehicle angular accelerations (or torques) are translated to gimbal rates 
is referred to as a steering law. 

V.2.2 CMG Steering Laws 

The equation of rigid-body angular motion for the vehicle (not including 
the CMG rotors and gimbals) is given by 

T v = lyffl + 0 I v u (4.2) 

where 

T v = 3 * 1 matrix of projections (onto the vehicle principal axes) of 
the net torque vector T v acting on the vehicle 

I = diagonal matrix of vehicle inertias about the principal axes 

u = 3 x 1 matrix of projections of the vehicle rate vector 

* 

u = 3 x 1 matrix of projections of the vehicle acceleration <j 
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a = 


-®z “y 


CO. 


-CO, 


X 


-co co 0 
y x 


C4.3) 


Similarly, the equation of motion for each CMG -rotor can be expressed in 
berms of angular momentum as 




+ 0H i 


(4.4) 


>rhere T p and H. are the 3 x 1 matrices of projections of the torque and angular 
K i 1 

nomentum vectors, respectively, for CMG number i. The net torque, T, produced 
by all six CMG's on the vehicle, neglecting gimbal inertias, is then given by 


T = - 2 T r 
i=l R i 


(4.5) 


ffhere 


= -H - QH 


H = 2 H 


(4.6) 


(4.7) 


(4.8) 


1=1 1 

By the law of conservation of angular momentum, 

+ ^ 

inhere 

H n = Vo + J T d dt ' 
o 

= net angular momentum of the vehicle plus the CMG 1 s 

X 

co = co at time t = o 
o 

= torque acting on the vehicle not caused by the CMG, 
i.e., T v = T + T d . 

The CMG system is normally initialized so that H is close to zero when u = 0; 
i.e., is normally close to zero. 

From the previous equations, the vehicle angular acceleration is given by 

« = -I T 1 [h + Q (Hn) + T d ] 


(4.9) 


= -I. 


H . 
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This approximation becomes an equality if co =0 and T, = 0. The -vehicle accel- 

* o » a 

eration, co, can thus be controlled by controlling H. If H can be controlled so 
that 


H = - 


Vc 


(4.10) 


■where co c is a commanded m, the torque produced on the vehicle by the CMG's is 
given by 


T = Lp _ + I (co - co 


V c 


(« “ a n) 


(4.11) 


- 1 „ 


where co n = E^, the vehicle rate not caused by the CMG system. For co s= 0, 

the CMG system thus inherently compensates for the centrifugal couple, 3L^o, 
produced when ca has components in two or more of the principal axes for vehicles 
with unequal inertias. 

By talcing the time-derivatives of the components of H, given in equations 
(4.1) for the 6-GAMS configuration, H can be expressed by 


H = A a 

where A is a 3 X 6 matrix of elements, a^, given by 


(4.12) 


a — h. S R C 
ij 2 P ctjL 


“a ■ h J ( s rj \ - C P C Vi \) 

hi ■ ( c r, \ + %\S) 


(4.13) 


for the stated configuration. Thus, the matrix A is a function of the gimbal 
angles a = ja-^, . “gj* 

* • 

H can therefore be controlled by controlling the gimbal rates, a. To 

control H according to equation (4.10), a solution of equation (4.12) is 
required. Since A is singular, any solution that exists is not unique. The 
specific solution that is adopted for a CMG system, even if it is only an ap- 
proximate solution, is commonly referred to as the CMG steering law. 

The selection of a steering law significantly affects the sizing (magni- 
tude of hj_) required for the CMG's in order for the system to be able to meet a 
specified angular momentum envelope requirement. The most prominent steering 


law is based on the pseudo -inverse 


2,3 


of A and is referred to as the 
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pseudo -inverse steering law. The significant feature of this inverse, denoted 
by A 1 ", is that it provides the unique minimum norm solution if a solution exists. 
In this case, A' can be computed by 



In minimizing the norm of the gimbal raxes, this steering law tends to 
emphasize the gimbal rates of the CMG’s that most efficiently provide momentum 
transfer in the commanded direction, thereby tending to maximize the net angular 
momentum envelope obtainable for the CMG configuration. Although other steering 
laws have been investigated (such as the superposition of solutions for three 
GMG's, and the transpose steering law which provides a very approximate solu- 
tion), they have been found to be either less efficient with no simplification 
in the computations, or very inefficient with considerable computational simpli- 
fication. It may be possible to improve on the pseudo -inverse steering law to 
increase the efficiency of angular momentum utilization even further, but a 
study to identify such a law Is beyond the scope of this project, Without fur- 
ther justification, the pseudo-inverse steering law is selected for the CMG con- 
trol computer described in this report. 

» 

The steering law computations thus provide gimbal rate commands, a Q , to 
the individual gimbal control systems in response to vehicle angular accelera- 

4 

tlon commands, co , In accordance with the following equation: 

a e = A+H = -A t I T d> c . (4.15) 

4.2.3 The Rate Gyro Configuration 

The vehicle rate information required for attitude control system stabil- 
ity and/or vehicle rate control is to be obtained from a set of six rate gyros 
configured to sense vehicle rates along the six dodecahedron axes shown in 
Figure 4-3. These axes correspond to the normals of the faces of a regular 
dodecahedron and have the unique property that the acuxe angles between any two 

axes are equal, given by 2&~63*4 degrees. This configuration is described by 
4 

Gilmore The significant feature of this sensor configuration is that the 

vehicle rates can be obtained from any three of the six dodecahedron rates, de- 
noted by r = ] r 1? ..., r^j. 
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The vehicle rates, <-», can he projected onto the dodecahedron axes to ob- 
tain the dodecahedron rates r "by 


r = E ce> 


(4.16) 


where 


E = 


r s 
-s 
c 
c 
0 
o 


0 

0 

s 

-s 

0 

c 


c 

c 

o 

0 

s 

-s 


(4.17) 


and 


s = sin b 


c = cos 6 




- 


1/2 


(4. IB) 


0.850 


If all six rate gyros are operating properly, the 'signals they generate, denoted 
by r s , are approximately equal to r. If a failure is detected in rate gyro num- 
ber i, however, this gyro is disengaged; i.e., r g ^ = 0. Let 


A = diag jAp 




(4.19) 


where' A^ = 1 if rate gyro number i is engaged, and A = 0 If it is disengaged. 
Then 

r s = Ar = Eco (4.20) 


where E = AE. To determine 'o from r 0 requires the solution of equation (4.20) 
Although S is a constant matrix for every state of A, there are b-2 states of A 
for which at least three of its elements are l’s; and there are. 20 states for 
which exactly three of its elements are l’s. Therefore, at least 20 different 
solutions are required in order to find co under all failure conditions (up to 
three failures). 
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Since any three rate gyros contain sufficient information to find co, only 
the 20 solutions are required to meet all failure conditions. By using all the 
rate gyros available, however, the effects of inaccuracies in the signals can be 
minimized. The pseudo -inverse solution given by equations (4.21) and (4.22) 
provides the solution for ao that minimizes the norm of Ar - r s (the "least- 
squares-fit" solution) 


where 

E 1- = (e T e) 1 E T 

= (e t ae) -1 (e 5 a) 



(4.21) 


(4.22) 


and where ffl s is the set of sensed and transformed signals for co. This solution 
is identical to a three-signals-only solution when all but the three respective 
elements of A are set to zero. 


The detection and isolation of failures in the rate gyro package is also 
a major task for the CMG control computer in order to generate A. Catastrophic 
failures can usually be detected by various types of monitors, but errors in the 
signals are not as easy to detect. Gilmore describes a set of parity equa- 

tions that isolate two gyro failures and detect a third failure. These equa- 
tions are presented in Table 4.1 under the notational convention adopted for 
this report. 
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TABLE 4.1 


Parity Equations for Rate Gyro Failure 
Detection and Isolation 
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Control Computer Description 


Hie principal function of the CMG control computer is to generate gimbal 
rate command signals to the six CMG's of the configuration described in psora- 
graph such that the vehicle responds in accordance with input command 

signals for vehicle rate and acceleration, using the rate signals provided by 
the six dodecahedron configured rate gyros described in paragraph V. 2.3« An ob- 
jective of the computer structure is to take advantage of the redundancy in both 
the CMG and rate gyro configurations to provide fail- operational performance, 
including failures m the computer. 

The control law for the vehicle rate control system, illustrated in 
Figure 4-4 can be written in the form 



•where 

<» cc _ vehicle acceleration commanded directly as an input 
to the rate control system 

» c " <» s = vehicle rate error 

= 3 x. 3 diagonal matrix of compensation functions 

The input signals :o cc and are provided either by a conxrol law for the ve- 
hicle attitude control system or by manual commands. The vehicle rate signal 
a> s is obtained by converting the six dodecahedron axes rate signals r g to 
vehicle axes signals. 

Since both the sensor signals r and the gimbal rate command signals 
a c are 6-tuples, the question has been raised as to whether there are any ad- 
vantages in performing the computations in six variables without first reducing 
the rate signals to the three-variable format. It has been suggested that by 
performing computations of data in a redundant format, perhaps some reliability 
advantage can be obtained with incremental computers, for which simultaneous 
computations are performed by separate components. The thought is that if one 
or more of such components should fail without affecting the remainder of the 
computations, the redundancy of the data and the computations would permit fail- 
operational performance. Investigations into this approach indicate that fail- 
operational performance is obtainable for certain component failures, but the 
complexity of such a system is considerably greater than for a system with re- 
dundancy at the overall computer level. In addition, there are many types of 
computer failures that do not permit fail-operational performance. Two complete 
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jomputers with failure monitoring or three such computers with majority voting 
provide greater reliability and appear to be less complex for reasons discussed 
subsequently. 

The vehicle accelerations in the dodecahedron axes for which rate gyro 
signals are available, Ar, can be related to gimbal rates by combining equations 
(4.9), (4.12) and (4.20) to obtain 


Ar = -Ha 


■there 

R = EI v _1 A 


(4.24) 


(4.25) 


md where the disturbance terms of equation (4.9) are neglected for simplicity. 
Che steering law in this case is required to solve equation (4.24) for a in 
terms of vehicle acceleration command signals r Q in the six-variate format, 
this command signal is given by 



(4.26) 


md for the vehicle rate control law of equation (4.22) 


r = E» + E G 3 (a. - co.) . (4.27) 

3y requiring that the three diagonal elements of G 3 be identical, thereby re- 
miring that the control system characteristic be identical in the three vehicle 
axes, the product E G J can be written as G E, where G is a diagonal 6x6 
natrix having diagonal elements identical to G 3 . Equation (4.27) can then be 
■mitten as 


c 



+ G 





(4.28) 


?igure 4-5 shows the structure for this control law. 


The steering law in this case is required to invert the 6x6 matrix R 
and the pseudo-inverse is selected for the reasons given in paragraph V. 2.2. 

Ohe simplest method for computing the pseudo -inverse of a 6 x 6 matrix of rank 3 
Ls to first factor it into the product of a 6 x 3 matrix M times a 3 x 6 
natrix If , which is always possible The pseudo-inverse S of R = MW can 

then be computed by 



(4.29) 
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(4.30) 


Some authors define the pseudo -inverse by equation (4.29). 7 Therefore 

S = N f M* 

S can be factored as follows to obtain 

s = (l T ' 1 A) t E 1 

= A t I v E 1 ' . (4.31) 


Although there are infinitely many ways of factoring R to obtain the unique 
S = R 7 , equation (4.30) indicates that the pseudo -inverse of this 6x6 matrix 
is equivalent; to a transformation of the six-variate acceleration commands f c 
to some three-axis format, followed by a transformation to the six-variable gim- 
bal rate commands a c< The question is whether implementation of the 6x6 
matrix S can provide any advantages m terms of simplicity or reliability over 
cascading the 6x3 and 3x6 matrix factors of S. 

Direct computation of the elements of S can provide limited fail- 
operational performance if such failures can be detected and isolated. By ex- 
panding the steering law computations 

a - -S r 
c c 

to obtain 

°cl = " S ll r cl ■ S 12 r c2 " S l6 r c6 

a c2 = _S 21 *cl “ S 22 *c2 “ •** " S 26 *c6 

• * • • 

a c6 = “ S 61 r cl “ S 62 r c2 " *** “ S 66 r c6 

it can be observed that the effects of a failure in the computation of the ele- 
ment, S^, can be nullified by letting A ^ = 0, thus causing r to be zero. 

The steering law is a function of A, and is therefore automatically modified to 
compensate for this step. Since the elements of column j in S are nullified by 
this procedure, isolated computation of the elements of this column is required, 
however, all 18 elements of A 7 must be computed for each of the six columns of 
S, necessitating six isolated computations of these elements, each element being 
n relatively complex function of the six gimbal angles. An increase in these 
computations, which constitute the bulk of the steering law computations, by a 


(4.32) 


(4.33) 
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factor of six represents a very significant increase m the complexity of the 
entire computer. Yet it provides fail-operational performance for only a 
limited set of failure types, namely those failures which can be isolated to the 
computation of the elements of a column. Further use of the rate gyro corres- 
ponding to the failed column is unfortunately also nullified by this procedure, 
thereby reducing the redundancy of the rate gyro package. 

To perform steering law computations on the six-variable signals also re- 
quires conversion of the three-variate input command signals, and co^, to r cc 

and r Q . In addition, six rate-loop compensation fillers are required instead of 
the three required with the three-variable system. 

• \ 

Based on these observations, it appears that parallel redundancy of the 
entire computer is preferred both from standpoints of complexity and reliability. 
The remainder of this report, therefore, considers only the three- variate struc- 
ture for the control law. 

Figure 4-6 shows the basic structure for the CMG control system m 
terms of its subcomputers, the CMG configuration, the vehicle, and the rate gyro 
package. Redundancy in the computation is not shown; this aspect is discussed 
in the following section. 

!+ * 2 • ? Fail-Operational Computation 

Fail-operational performance for a computation function can be obtained 
with two or more complete channels of computation plus monitors that are able 
to detect and identify a channel failure. In some cases, uhe only reliable 
technique for identifying a failure is to "vote" between at least three channels 
for a single failure, at least four channels for two failures, and so on. If 
it is possible to reliably detect a failure in a single channel by a simple 
method, one less channel of computation is required in comparison with the 
voting method. It is therefore worthwhile to investigate techniques for failure 
monitoring. 

The steering law computation and the dodecahedron inversion both consist 
of the inversion of matrices that are very much simpler to compute in the for- 
ward direction than to invert. By performing the forward transformation M on 
the output of an inversion transformation M* and comparing the resulting out- 
put with the original input, an indication of a failure m either of these 
transformations is obtained. This principle is illustrated in Figure 4-7. 

The scope of this study does not permit comparisons of various methods 
for detecting and isolating failures, and the approach suggested here requires 
further study. 
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k.3 STEERING LAW COMPUTATION 

Since the steering law computation is by far the most complex function to 
be performed by the CMG control computer, most of the effort under this program 
has been devoted to this function. Paragraph *1.3.1 presents the equations to be 
computed by the steering law computer, paragraph *1.3. 2 describes the principles 
of incremental computation, and the remaining paragraphs describe how the equa- 
tions can be. implemented with incremental computation. 

*1.3.1 Equations and Computer Structure 

The steering law computer is required to generate glmbal rate command 
signals a c to the individual gimbal control systems of the CMG* s in response 
to vehicle angular acceleration signals u e in accordance with equations given 
in paragraph * 1 . 2 . 2 and repeated below for convenience 



I . I . and I are the inertias of the vehicle about its principal axes, and 
x 5 y’ z 

[a* = A~ (AA 1 )^ 1 ] is the pseudo -inverse of A. 
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To simplify notations In the subsequent development, the following defi- 
nitions are made. 


B = 

AA' T 

(4.34) 

C = 

ad 3 B 

(4.35) 

II 

O 

K) 

det B 

(4.36) 

D = 

a t c 

(4.37) 

u = 

(%)M ■ 

(4.38) 

A* 

= D (t) 

(4.39) 


and 


d c = -DU 


(4.40) 


By leaving A" 5 ” in the factored form, as in equation (4.39) only three divisions 
by the scalar d are reaulred comnared to 18 such divisions if the elements of 
were to be computed. Since B is only a 3x3 matrix, the adjoint method for 
inverting B is used. The steering law computer can now be structured of sub- 
computers related to the above defined variables so that they can be discussed 
separately. Figure 4-8 illustrates this structure of the steering law 
computer. 


The A computer computes the elements of A given by 

a. = h, S C 
ij 3 P a.j 

= h , (s S - C c C ) 

2j 3 \ a j e Tj ) 

a,. = -h, (c S + C_ S C \ . 

33 i [ r 3 a 3 P Vj a 3 ) 


(4.41) 
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( 4 . 42 ) 




Element 


may be factored into 


a i3 “ p iJ ff J 


where 

n . . = (m. . S + n . - C \ 

Pij { ij a d ±3 a 3 ) 


( 4 . 43 ) 


and. where is the speed of rotor j in revolutions per second. The constants 
and n^. are given by 


il 

sf 

0 

n lj ~ 2nI R% 


m 2j - 

2 n Ix>S 
R y } 

n 2j = “ 2TTl R C 3 C y 3 

(4.44) 

n 

■n 

a" 1 

-2 n I R C 
R y 3 

n 3j = - 27cI R c p s r 



where is the inertia of the CMG rotor about its spin axis. For the 6-GAMS 
configuration described in paragraph 4.2.1, the y angles are 0, 60, 120, 180, 
240, and 300 degrees so js j and jC j are specified, but the gimbal axis uilt 

angle g depends on the requirements for the net angular momentum envelope. 


Thus, the A computer for each of the 18 elements ] a ij | performs the 

scaling and addition of the input signals, S and C , as indicated by ectuation 

a J a j 

(4.3-10) and multiplies the result by the input signal for rotor speed, a . 

J 


To help visualize the complexity associated with each of the subcomputers 
of the steering law computer, the related matrix expressions are expanded in the 
following discussion. Since B = AA + is a symmetric 3x3 matrix, only six elements 
such as those making up the upper triangle of B, are required. These elements 
are 


“'ll 


J 12 


13 


22 


23 


33 


2 2 2 2 2 2 
a ll + a l2 +a i3 + a i4 + a l5 + a l6 

a ll a 21 + a 12 a 22 + a .13 a 23 + a l4 a 24 + a 15 a 25 + a l6 a 26 

a ll a 31 + a l2 a 32 + a 13 a 33 + a l4 a 34 + a l5 a 35 + a l6 a 36 
a 2i 2 + a 22 2 + a 23 2 + a 2lf 2 + a 25 2 + a 26 2 

a 21 a 31 + a 22 a 32 + a 23 a 33 + a 24 a 34 + a 25 a 35 + a 26 a 36 

2 2 2 2 2 2 
a 31 + a 32 + a 33 + a 34 + a 35 + a 36 


(4.45) 
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The C computer evaluates the upper triangle elements of C = adj B since 
this matrix is also symetric = 

C 11 = b 22 b 33 “ b 23 2 

c 12 = b l3 b 23 " b l2 b 33 


c 13 = b 12 b 23 " b 13 b 22 


c 22 ' b ii b 33 “ b l3* 


c 23 = b 12 b 13 " b ll b 23 


c 33 = b ll b 22 " b 12 ‘ 


Elements of B and C are used to compute d Q : 


( 4 . 46 ) 


- aX 


d o = 

b ll c ll + 

b 12 c 12 + 

b 13°13 

c are 

computed 

as follows : 

d li = 

a ll c ll + 

a 21 c 12 + 

a 31 c 13 

fi 12 = 

a ll c 12 + 

a 21 e 22 + 

a 3l c 23 

d 13 = 

a ll C 13 + 

a 21 c 23 + 

a 31 c 33 

d 21 = 

a 12 c ll + 

a 22 c 12 + 

a 32 c l3 

d 22 = 

a 12 C 12 + 

a 22 c 22 + 

a 32 c 23 

d 23 = 

a 12 c 13 + 

a 22 c 23 + 

a 32 c 33 

d 31 = 

a 13 c ll + 

a 23 C 12 + 

a 33 c l3 

d 32 = 

a l3 c 12 + 

a 23 C 22 + 

a 33 c 23 

d 33 = 

a 13 c 13 + 

a 23 c 23 + 

a 33 c 33 

d 4l “ 

a l4 c ll + 

a 24 c 12 + 

a 3 1 + c 13 

d 42 “ 

a l4 c 12 + 

a 24 c 22 + 

a 34 c 23 

\ 3 = 

a l4 c 13 + 

a 24 c 23 + 

a 34 c 33 


( 4 . 47 ) 


( 4 . 48 ) 
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a 51 = a 15 c H + a 25 c 12 + a 35 c l3 
a ?2 = a 15 C 12 + a 25 C 22 + a 35 C 23 
a 53 = a 15 c 13 + a 25 c 23 + a 35 c 33 
61 a l6 c ll + a 2 6 c 12 + a 36 c l3 
a 62 = a 16 c 12 + a 26 c 22 + a 36 c 23 
d 63 = a l6 c 13 + a 26 c 23 + a 36 e 33 * 


The only divisions required in the steering law computer are in the u 
jomputer. The elements of this 3x1 matrix are 

■^cx 

u i = ~aT~ 


u _ h&x. 

w-o x 


( 4 . 49 ) 




u 3 ~ <L 


The a c computer then provides the outputs of the steering law computer by 
performing the following computations' 

“■el = a il a l + fl 12 u 2 + a l3 u 3 


d c2 " a 21 u l + a 22 u 2 + d 23 U 3 


“c3 “ d 31 u l + a 32 u 2 + a 33 u 3 


( 4 . 50 ) 


“cH = a 4l u l + <V. U 2 + a ^3 u 3 


“c? ~ d 51 u l ^ d 52 u 2 + a 53 u 3 
“c6 = a 61 a l + d 62 u 2 + d 63 u 3 • 
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A summary of the types and numbers of mathematical operations required 
for the subcomputations of the steering law computer is presented in Table 4.2. 


TABLE 4.2 

Mathematic Operations for the 
Steering Law Computer 


Subcomputer 

Additions 

Mul - tipi i cations 

Divisions 

A 

18 

18 

0 

B 

30 

36 

0 

C 

6 

12 

0 

d o 

2 

3 

0 

D 

36 

54 

0 

u 

0 

0 

3 

& 

c 

12 

>N 

18 

0 

Total 

104 

l4l 

3 
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k 3-2 Incremental Commutation 


An incremental computer is a special purpose digital computer that has a 
lumber of features which make it attractive for solving equations of the type 
required for the p seudo -inver se steering law. In contrast with the general 
purpose type digital computer, which performs entire computations during each 
computation cycle, an incremental computer updates previous computations. Such 
computations are generally much simpler and are performed simultaneously on the 
lumerous variables of the problem. The time required to complete each computa- 
tion cycle is therefore much less than for a general purpose computer, permit- 
ting higher sampling rates. 

In an incremental computer, the variables of the computations are stored 
Ln binary registers, referred to as Y registers. The. content (numerical value) 
stored in a Y register is increased by one when so commanded by a pulse, AY, 
representing an incremental increase in the variable Y. The AY pulse can also 
be negative, in which case the Y register is decreased by one. In addition, the 
content of a Y register is added to the content of an R register when so com- 
nanded by a positive nX pulse, or it is subtracted from the R register for a 
negative AX pulse. 

Figure ^.3-2 shows the schematic symbols adopted for the incremental com- 
puter elements. Since the R register has the same capacity as the Y register, 
repeated AX pulses of the same sign will cause the R register to overflow. Each 
time the R register overflows in the positive direction, a positive AZ pulse is 
generated; if the overflow Is In the negative direction, a negative AZ pulse is 
generated. This AZ pulse can then serve as a aX pulse or a AY pulse for other Y 
and R registers. Let 


Z(n) = S AZ(1) (4.51) 

1=1 


where AZ(i) is the AZ pulse produced by the i'th AX pulse, AX(i) [AZ(i) may he 
0, +1, or -l] t£ Let Y(i) and R(I) be the contents of the Y and R registers, re- 
spectively, at the- occurrence of AX(i), and let c be the capacity of the regis- 
ters. Then 


cZ(n) + R (n) = R(l) + S Y(i) AX(i) . (4.52) 

i=l 


Therefore, for Y(i) > > 1, 


Z(n) 


i n 

- 2 Y(l) AX(I) • 
c 1=1 


(4.53) 
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For the case where Y(I) is a function, of X(i), 


ZCn) ~ 


1 

c 


X(n) 

Y[X(i)] dX(i) . 


(4.54) 


This computation structure may thus he used to integrate functions of indepen- 
dent variables, and combinations of such elements may be used to generate func- 
tions that are solutions of differential equations. Computers structured by 
interconnecting integrators of this type are commonly referred to as DDA’s 
(digital d Iff erential analyzers). 


The computations for the steering law require only additions, multiplica- 
tions, and divisions. The sum of several variables can be obtained by feeding 
their respective A pulses to a common Y register and properly controlling their 
pulse times so that they do not occur simultaneously. Such timing control can 
be implemented by various methods, and will be discussed later. 

Multiplication is readily obtained by noting that 

a (V 2 ) = V T 2 4 Vh (4.55) 

Therefore, let AX^ = AY 2 and AXg = AYp and let the two Y registers add into a 
common R register as shown in Figure 4-10. Then 

AZ ~ c (h AY 2 4 Y 2 M l) = i A ( Y 1 Y 2 ) • (4.55) 

The accuracy of this computation improves with the size c of the registers. 

It can also be improved oy proper sequencing of the operations when AY^ and AY 2 
occur simultaneously. Investigations of such methods are beyond the scope of 
this study. 

4.3.3 AA Computation 

The AA computer generates increments jAa^j to increase or decrease the 
values in the set of 18 registers that store the elements of the A matrix, given 
by equations (4.42), (4.43) and (4.44). Column j of A corresponds to CMG 

number j. The incremental computer structure for the three elements of this 
column is shown in Figure 4-11. 


The Y register corresponding to some variable, y, is identified by this 
variable, and It stores the numerical value 


„4.^ ,™4 4V44. K y y where K y Is a scale factor. 

The R register that generates Ay is Identified by y. Constant numbers are 
represented by the oval-shaped areas in Figure 4-11. For example, when a posi- 
tive pulse for Ac ^ occurs, the Y register containing K Q is increased by the 

constant number K . 

3 


'3 
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FIG. 4-9 Basic incremental computer elements. 


r 1 ^ 



FIG. 4-10 Incremental multiplier. 
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The rotor speed o ^ (revolutions per second) can be held constant by the 
spin motor electronics to some degree of accuracy, but high accuracy may be dif- 
ficult to achieve. If the rotor speed is sensed and Included in the steering 
law computations, rotor speed control is not required at all, permitting simple 
spin motor electronics. The rotor speed is therefore included in the steering 
law computations. 


The rotor speed is computed by counting revolutions 0^ over a period of 
time Tg.. A pulse A0^ produced for each rotor revolution increases the 
counter by one. A pulse AT a . is produced each T a seconds by some clock. 
must be quite long, such as several seconds, to permit sufficient resolution in 
ffj. When the AT^ pulse occurs, the 0^ register is first set to Kg 0 j, where 0 j 

is the present number in the counter and where Kg = K 0 /T a (chosen to be a 
power of 2 for simplicity). The counter is then reset to zero. Next, Kg 0, and 
K fl d, are compared by adding the first to, and subtracting the second from, the 

e register. If this difference is positive, a positive ac, pulse is generated 

J J 

and causes the a, register to increase by K , or vice versa for a negative 

J 1 

aulse. 0 


To obtain maximum accuracy in incremental computations, the Y registers 

should be scaled so low that they are as full as possible without overflowing. 

Therefore, K_ is selected so that K„ times maximum o , is close to but does not 

"3 “j 3 

exceed the register capacity c. 

T a |or Kg = K 0 /T should be selected to minimize the error in cr^. If 

T ff is very large, the maximum error in due to resolution, given by 

or j/0j = 1/T^, is small. However, the maximum error due to a change in ffj during 

the T„ interval. cr_T„ (where = max o, assumed constant in this interval), is 

o J in o m 7 — — 

large. The T a for which these errors are equal is given by T ff = 1/y • 


Incremental signals for sin and cos Oj can be obtained by numerous 
methods, but at the present state of the art a set of six resolvers with a 
single time-shared A/D converter for all six gimbal angles appears to be the 
test candidate. Selection of a preferred method Is not performed under this 
study. The scale factor for the registers in the A/D converter that store S a 


Is K_; i.e., a AS pulse corresponds to a change of 1/K 


sc J 


sc in S a^ 


etc. 


3 


and 
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When a positive AS a pulse or AC a pulse occurs, the registers are 
increased by the respective constants, m^ and n^ . If this produces a positive 
overflow of a p^. register, the resulting Ap^ pulse causes p^ to increase by 
one, and cr^. to be added to a^ . Overflows in the a^ registers then produce the 

output pulses of the AA computer, which are held in flip-flops until called for 
in the next computation cycle. Additions to the registers are also produced 

by a Ao^ pulse. 


Although there are many ways to perform the functional operations symbol- 
ized in Figure 4-11, considerable savings in hardware may be realized by 
serializing many of the variables on a single, long, circulating register. 

Table 4.3 describes one method for sequencing the additions into the a^ 

registers. All 18 variables | Pi j f are ser f a H ze( I on a single, long register 

in the sequence shown in the table, and the size variables are serialized 

on a single register such that exactly three cycles of the 0 register coincide 
with one cycle of the p register, and such that the least significant bits of 
°1’ Pll’ an ^ E ij coincide, etc, as shown. The 18 a registers are the same 

length as each word of p and a. During word interval 1, the p register adds 
into the a^ register if = 1 (or subtracts if A0-j_ = ~1 )j the c register adds 

into a^g if Ap_^g = 1 (or subtracts if A P^g = - 1)3 etc, as shown in the table. 

It is also possible to serialize the a., variable with increased logic com- 
plexity and longer cycle time. 
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TABLE 4 . 3 

S*t 

A Computation Sequence 


Interval: 

1 

2 

“ ~ 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

l4 

15 

16 

17 

18 

pi 

Pn 

P 21 

p 3l 

Pl2 

p 22 

p 32 

p 13 

P 23 

p 33 

Pl4 

P 24 

P34 

p l5 

p 25 

p 25 

p l6 

p 26 

p 36 

<Jt 

*i 

°2 

ff 3 

“4 

°5 

a 6 

c l 

*2 

°3 

°4 

*5 

°6 


°2 

°3 

°4 

a 5 

tf 6 

"n 

a 21 

“31 

*12 

a 22 

*32 

*13 

a 23 

a 33 

*l4 

*24 

*34 

*15 

a 25 

*35 

*16 

*26 

*36 

pia 1 

aa Pl6 
**P 26 
ffAP 36 

P A °1 

P*°l 

P A< *2 

P a0 2 

P A0 2 

CTA Pll 

ffAp 21 

PAp 3 i 

P^ 0 3 

CT * p 12 
**P 22 

dip 32 
P AC 3 

OAP13 

<JAp 23 

P a °3 

p A c, 4 

P a °4 

pA CTlf 

p a<j 5 



\ 

fl Ao 5 

“ Ap 23 
°* p 33 

p A o 5 



CA Pl4 

°*P24 

°* p 34 

P ao 6 

0A P 1? 

«Ap 25 

ffAp 35 

pio 6 

pA<J 6 



The Aa^ 

aB computer). 


pulses accumulate in the register (which is part of the 
The number stored in this register is given by [per equation 


(4.56)3 

a. (ho f (Ko.) . 

13 c V P'iD) V a 3 / 


(4.57) 


Therefore, since = p.^ 

quired to satisfy 


[equation 74.42)], 




the scale factors are re- 


(4.58) 


In addition, 


the Y registers cannot 


K 




'13 


K 


13 


be allowed to overflow, 
max jpij| c 
max £ c 

max t a ijl £ G 


Therefore, 


(4.59) 

(4.60) 

(4.61) 


If more than one Aa^ pulse is permitted to occur for each computation 
cycle, the complexity of the Aa^ adder must be increased significantly. To 
avoid this situation, let 

K max jp.,,j + 3C max |cr. I s c . (4.62) 

p i3 i 13 l 1 

Maximum accuracy is obtained with the scale factors as large as possible, sub- 
ject to the above constraints. 


Let the constants m.y and have the same scale factor, K ffl 


ij 


_ 1 




ij 13 

£ 3 




sc 


/ m. s + n -c \ 

( iJa 3 1Ja jj 


13 


Then 


(4.63) 


Since p, . = m. .s + n. .c [equation (4.43)], the scale factors must satisfy 
M ij ij 13 a 1 


K = K m, fsc 

P U — if— 


(4.64) 


Also, to prevent more than one Ap^ pulse per computation cycle, 

(i"i 3 r hii) s 0 


(4.65) 
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4.3.4 AB Computation 

The AB computer calculates the increments for the six upper triangle 
elements of the B matrix given in equations (4.45) . Figure 4-12 shows the 
structure of this computer, which has 18 Y registers to store the A matrix 
elements and six R registers to generate the Ab^ pulses. 

Each of the three diagonal elements of B is the sum of the squares of 
the six elements in a row of A. By letting aX = AY (see Figure 4-9), we 
obtain 


AZ 



fe My 2 ) . 


(4.66) 


Therefore, by adding K a . . to K. * b. . when Aa = +1 for each j, the 

a i3 13 D ii 11 *3 

content of the Y register that accumulates the Ab^ pulses will be 



(4.67) 


In order for 


b ii - a lj 


the scale factors for the diagonal elements of B must satisfy 

b ii 2c 


(4.6B) 


(4.69) 


where K = K for all j. 
a i a i3 


For the off-diagonal elements, multiplication is performed as described 
previously. In this case, 

6 

X E 




c 


K=1 


1,,K 9- 

ik a^ jk 


(4.70) 


Therefore, 




c 


i * 1 


(4.71) 
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The AC computer, described in paragraph ^.3.?, requires that all ele- 
ments of B have the same scale factor. Therefore, let IC = K. for all i and 

D 

j* The result is that must equal X^/2, which may be accomplished either 

3_i 

by adding twice to the registers for each Aa^ or by delaying the addi- 
tions by one bit, thus doubling the number added to the register. 

In order to ensure that, at most, one Ab.^ pulse is produced during each 

computation cycle, the following constraint is placed on K : 

« 

6 

K a K=1 j a ik| + max | a jk( * c for a11 j ’ C4<72) 


This constraint limits to a lower value than the constraint given by equa- 
tion (^.11). Since K = K for all i and j and since it is reasonable to 

a ij a 

let - K ff for all j (all wheels have the same speed), equation (4.50) in- 
*3 

dicates that K = K for all i and j. Equation (4.58) is therefore 
p ij p 
replaced with - 


K 


K K 
P p 


(4.73) 


Table 4.4 shows a computation sequence for the additions to the b^j 
registers where the variables are stored serially on a single circulating 

register in the order shown. Each register recycles in the period of a 
one-word interval of the A register. During interval 1, twice the number in 
the A register during that interval is added to the register if Aa xi = 1. 
The A register also adds to the register if &a 2 j_ ~ 1 during this interval, 

and to the b^ register if Aa^ = 1, All the computations of the AB computer 
are completed in one cycle of the A register. 
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The AC computer, illustrated in Figure 4.13, computes increments for the 
adjoint of the B matrix* given by equations (4.46). As in the AB computer, 
when squaring computation is performed, the addition to an E register is doubled. 
This can be accomplished by delaying the number to be added by one bit time. 

Ehe doubled addition is signified in the figure by a "2 11 placed next to the AX 
arrowhead which indicates addition of X to E. As fox' the previous computers, 


the scale factors must satisfy 



(4.74) 


By requiring that 

2 maxjb^l £ c (4.75) 

Where the summation is taken over the four terms b^ of each of equations 
(4.46), only one Ac pulse will be produced during each computation cycle. 

Table 4.5 presents a computation sequence for the additions to the 
c registers for this computation, where the elements of B are stored serially on 
a single circulating register in the sequence shown. 


TABLE 4.5 

G and d Q Computation Sequences 
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,V.3*6 Ad Q Computation 


Computation of the increments for d Q , the determinant at the B matrix 
[equation (4.47)] is performed by adding elements of B and C to a single 
E register, as shown in Figure 4-14. The sequence of these additions is 
shown in Table 4.5, where' the elements of C are stored on a single register 
in the sequence shown. The scale factor for d Q is given by 



(4.76) 


and by constraining and K Q so that 

3 . 

K^K q s ntaxjb |+ maxjc-^f £ c (4.77) 

u 1 

with only one Ad c pulse produced during each computation cycle. 


lf.3.7 AD Computation 


Figure 4-15 Illustrates computation of the Increments for the 18 D 
elements given by equations (4 48) . Similarly to the previous computers, the 
scale factors for the elements of D must satisfy 


K, 


K K 
a c 


‘1J 


= K, 


(4.78) 


Therefore, they must all be equal. By requiring that 

3 

K a K c = S max|a ki |+ max| k ^ 1 s c for all l, j (4.79) 

1c — 1 

only one Ad pulse is produced during each computation cycle. 

/v 

Table 4.6 presents a computation sequence for additions to the 'd 
registers where the elements of A and C are stored serially on respective single 
registers in the relative sequence shown. 
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Ac,, AC* ' C j; AC 



FIG. 4-15 AD computer. 
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TABLE 4.6 

D Computation Sequence 


Intermix 

X 

2 

3 


5 

6 

7 

8 

9 

10 

11 

12 

13 

lh 

15 

16 

1? 

18 

At 

*11 

“21 

*31 

*12 

*22 

"32 

*13 

*23 

'33 

*l*f 

‘24 

'34 

‘15 

'25 

*35 

*16 

'26 

*36 

Ci 

C 11 

'12 

e 13 

'22 

'23 

°33 

'n 

'12 

'13 

°22 

'23 

'33 

El 

'12 

C 13 

c 22 

'23 

c 33 

*n 

“'ll 

“'12 

“'13 






““31 










*12 

“'12 

“*22 

“'23 







“‘21 

“ 2 31 








*13 

“'13 

“'23 

“'33 






“'ll 


“*21 

“‘31 







*21 




“'11 

“'12 

“'13 

“*12 

“'22 

“'32 










*22 




“'12 

“'22 

“°23 


“'l2 


“'22 

“*32 








*23 




“'13 

“'23 

“'33 



“‘12 


“'22 

“'32 







*31 







“'ll 

“'12 

“'13 




ca» 13 

“'23 

“'33 




*32 







“'12 

“'22 

“ c 23 





“'13 


“'23 

“*33 


S 33 







“'13 

“'23 

“'33 






“‘13 


C1 ‘23 

“'33 

*41 










“'ll 

“'12 

“'13 

“‘14 

“*24 





*42 










“ C 12 

AAc 22 

A.c 23 


“*14 


“*24 

* *‘34 


*43 










“'13 

“'23 

“'33 



“*14 


“*24 

“*34 

*51 

CAa l5 

“ a 25 

“ a 35 










“'ll 

“'12 

“'13 




*52 





“ a 35 








“ c 12 

“'22 

“'23 




*53 



“*15 


“*25 

“'35 







“'13 

“ c 23 

“'33 




*61 

“'l6 

“ a 26 

CAa 36 













“'ll 

“'12 

“'13 

*62 


“ a 16 


“*26 

“'36 











“'12 

“'22 

“'23 

*63 



C1 *16 


“ a 26 

“*36 










“'13 

“'23 

“'33 
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4-- 3 • 8 au Computation 

The An computer performs the divisions of equations (4.49) by multipli- 
cations in the feedback paths, as illustrated in Figure 4-16. 

When the number stored in the register is positive, the element 
identified by xhe letter S in the figure (and inappropriately but convention- 
ally termed a "servo")) produces positive Au pulses at the computation cycle 
rate. When is negative, the servo produces negative Au. pulses. No pulses 
are produced when is zero. 

When a positive Aco cx pulse occurs, the number Kj I x , which is scaled 

to be nearly as large as c, is added to the £-^ register. Similarly, a positive 
Ad pulse causes K u, to he added to s', . The resulting Au, pulses produced 

O 11 1 -i- ^ 

cause d Q to be repeatedly subtracted from s-^ until is aero. (The register 
o 

may not go to zero exactly, in which case the servo will generate alternate pos- 
itive and negative pulses.) 

The change in due to pulses or Ad Q pulses can be expressed as 


As l = K I J x A “cx " Vl id o ' K d d o Au l 


= A 


(v*) 


K* co 
“x cx 


(Vi) (v») 


Since 0, and with proper initialization of the d Q and u^ registers, 


(4.80) 

(4.81) 


U, 


K t K, 
J X -°x 


X x j:> cx 


(4.82) 


To satisfy equations (4.49), the scale factors must therefore satisfy 


£ X - I 

*u K d 

o 


(4.83) 


Similar requirements apply to the y and z axes. 

To ensure that the e-^ register does not overflow, the scale factors 
Should also satisfy 


K I X x + ^ max U 1 + K d maX d c 


(4.84) 
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FIG. 4-16 Au computer. 
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Table 4.7 presents a computation sequence for e 2 , and where u-^, 
Ugj and u^ are stored serially on a single register as shown, and where d Q is 
stored on a single register that recycles at the rate of each element of u. The 
computation is completed in two cycles of the u register. 


TABLE 4.7 

s Computation Sequence 


Interval: 

1 

2 

3 



6 

u: 

U 1 

u 2 

u 3 

U 1 

u 2 

u 3 

a : 

d„ 

d. 

d 

d 

d 

d 

0 

O 

0 

0 

0 

0 

0 


uAd 

0 



d o AUi 



/V 


UAd 



d Au„ 


2 


0 



0 2 


S 3' 



UAd Q 



d Q Au 3 


J +. 3* 9 Aa c Computation 

The Aa c computer, illustrated in Figure 4 - 17 , computes increments for 
the desired outputs of the steering law computer [equations (4.50)]. As in 
the previous cases, the scale factor for the registers that accumulate the 
Ac c pulses (not shown) must satisfy 


K = 


¥u 

c 


(4.85) 


By requiring that 


K 2 max 

a . -r. 

0=1 


d ij 


max 


U-! 


< c for each 1 


(4.86) 


only one a Q register overflow can occur during each computation cycle. 


Table 4.8 shows a computation sequence for the a c registers, where the 
elements of D and u are each stored serially on single long registers in the 
sequence shown. 


Since u is most likely desired in analog form for the CMG gimhal control 
system, the pulses may "be fed Into counters that are parts of a D/A conver- 
sion system. 
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TABLE 4.0 


a c Computat ion Sequence 


Intervals 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

1 ? 

16 

17 

18 

Di 

d ll 

d 12 

d 13 

d 21 

d 22 

d 2 3 

d 31 

d 32 

d 33 

d 4l 

d 42 

d 4 3 

d 5i 

d 52 

d 53 

d 6 l 

d 62 

d 63 

u: 

U 1 

u 2 

u 3 

U 1 

u 2 

u 3 

U 1 

u 2 

u 3 

U 1 

u 2 

u 3 

u i 

u 2 

u 3 

U 1 

u 2 

u 3 

“cl 

Dau 2 

DAu 2 

DAu 3 

uad ll 

Ufid 12 

UAd 13 













a c 2 


uAd 22 


DAu^ 

Dau 2 

DAU 3 

u * d 21 












■«3 



u4d 33 

uAd 3 i 

uAd^g 


DAu 1 

DAu 2 

D au 3 










“c4 




uAd^ 

u&D 42 

UAd^ 




DAu^ 

dau 2 

DAu 3 







“c5 





md !fe 

UAd ?3 

UAd ?1 






DA^ 

DAu 2 

DAu 3 




“c 6 







eAd 63 

UAd 6l 

uAd 62 








DAu l 

Dau 2 

DAu 3 



4.4 OTHER COMPUTATIONS 


4.4.1 Triad Conversion 

The triad converter changes the six rate signals r g = j r s i>***5 r s6 [ 
from the dodecahedron configured rate gyros, described in paragraph 4.2.3, to 
vehicle rates or = jaw, a)„ v , co_ _ J . The equations to be computed are given by 
equations (4.21) and (4.22) and are repeated here for convenience : 


where 

E + = (E^E)" 1 E T 


o c 
o c 
s o 
-s o 
c s 
c -s 

A = diag |a^, •**> 

s = sin b « 0.526 
C = COS 6 sk 0.850 



If rate gyro number 1 is used, A^ = 1; if it is switched off because it 
has failed, A^ = 0. Let 


F = E T E 

(4.87) 

P - adj F 

(4.08) 

Q = PE T 

(4.89) 

q Q = det F . 

(4.90) 

oi — ~ — Qr„ 

(4.91) 

3 % 5 
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Since E changes with A, the Q matrix can he recomputed by performing these 
computations each time A changes, 'this method requires minimum memory and is 
probably the preferred method with a general purpose computer. The elements of 
Q can also be expressed in terms of the elements of A. This method is pre- 
ferred for incremental and analog computation, and is described as follows. 


*1 = 

„ 4 
as c 

( 4 . 92 ) 

n 2 = 

2 s^ + (i-j_ 

( 4 . 93 ) 

^3 = 

»*1 + ^2 

( 4 . 94 ) 

^4 = 

*2 + ^3 

( 4 . 95 ) 

A i jk 

= A i A j A k • 

( 4 . 96 ) 


Then the 18 elements of Q are given by the following expressions: 


q ll 5 

= !-2 A 123 H 

h ^2 A 124 4 

! 1 4 A i25 


+11 4 A 126 

+ tx 2 A 13? 

“ p 1 a 136 


"^145 

+ ^2 A 146 

+ 2 ^ 2 A 156 

q 12 5 

" _Ii 2 A 123 

- ^2 A 124 

- ^^125 


-tl 4 A l26 

+ ^1 A 235 

- ^ 2 A 236 


_tl 2 A 245 

+ ^1 A 246 

+ 2 ^ 2 A 256 


q 13 ~ ^3 A 134 + ^4 a 135 + li 3 A 136 
+| 1 3 A 234 + ^3 A 235 + ^236 
+ 11 1 A 346 + ^1 A 346 + 2 ^3 A 356 
q l4 = ^3 a 134 + ^3 A l45 + ^4 A l46 
+ > l 3 A 234 + ^4 A 24? + ^3 A 246 
+ ' i l A 345 + ^1 A 346 + 2 ^3 A 456 


- 82 - 



q i$ = -,1 3 A 135 + ^145 “ ^3 A i56 
-^ 2 A 235 + ^ 3^45 + ^ 3 A 256 
~^2 A 356 + ^2 A 456 
q 16 = _fl 2 A 136 + ^3 A l46 + ^3 A 156 
“**3 A 236 + Ji 2 A 2 ! t6 “ ^ 3 A 2?6 
_ ^2 A 356 + (1 2 A 456 
q 21 = _ll 3 A 123 + ^3 A 12lf " lt 2 A 125 
+li 2 A 126 ~ fl 3 A 135 + ^2 A 135 
-tl 2 A lif5 + ^3 A l46 

q 22 = ^3^123 “ ^3 A 124 " |l 2 A 12? (4.97) 

+(i 2 A 126 ~ ^2 A 235 + ^3 A 236 
-M- 3 a 2 45 + ^2 A 246 
q 23 = 2fl 2 A 123 + ^134 + ,1 2 A 135 
_,1 1 A 136 + ^4 A 234 " ^l A 235 
+ ^2 A 236 + !1 2 A 345 + ^2 A 346 
q 24 = “ 2fl 2 A 124 " **4*134 " * 1 1 a i 45 
_ti 2 A l46 “ **4 A 234 " fi 2 A 24? 

+lt l A 246 " ^2 A 345 “ J1 2 A 346 
q 2j = 2 **3 A 125 + * 1 4 a 135 + p, 3 A l4? 

+|i i A l56 + |i 3 A 235 + F 4 A 24? 

+ **1 a 256 + ^3 A 356 + !i 3 A 4J6 
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q 26 = 2 *3 A 126 + *3 A 136 + *4 A l46 
+ ^3 A 156 + ^4 A 236 + *3 A 246 
+ti l A 2?6 + *3 A 3?6 + ^3 A 456 
q 3l = |i l A l23 + ti l A 124 + ^3 A 135 
+ *3 A 126 + 2|1 3 A 134 + ^135 
“ ll 3 A 136 + |l 3 A l45' + *4 A l46 
q 32 = *1 A 123 + [i l A 124 + |a 3 A 125 

+ll 3 A 126 + 2 ^3 A 234 + • l 3 A 235 
+ *4 A 236 + *4 A 245 + 11 3 A 246 

C4.97) 

q 33 = " tl '2 A 134 " fl 3 A 135 “ ft 2 A 136 
+ *2 A 234 + |1 2 A 23? + ll 3 A 236 

-^3 A 3 i +5 + 

q 34 = ~*2 A 134 " *2 A l45 " ii 3 A l46 
+ *2 A 234 + ^3 A 245 + ^2 A 246 
+|i 3 A 3 l !-5 “ ^3 A 346 
q 35 = *2 A 135 ‘ ^1 A 14? + *2 A 156 
“* i l A 235 + *2 A 245 + *2 A 2?6 
+2^2 a 3 4 ? + *4*356 + * t 4 A 456 
q 36 = **1*136 " **2 A l46 “ u 2 A 1J6 

”*2 A 236 + |1 1 A 246 “ *2 A 256 
-2ji 2 A346 - *4*356 - *4*456 
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( 4 . 98 ) 


To express q in a similar maimer, let 

l; = 4»V = (.3 - s-’) 2 

„ 6 - Vs 2 = 4 - (o3 + ,3)2 . 

Then 

q o = jl 5 A 123 + >"5 A 12V + u 6 A 125 
+ V - 6 A 126 + 11 6*134 + ^6 A 135 
+ ^5 A 136 + **5*14? + ( 1 6 A l46 
+ > i 5 A l56 + ^6 a 234 + ^5 A 2 35 
+ ^ 6 A 236 + ^ 6 A 245 + 1 * 5*246 
+ ^5 a 2?6 + il 5 A 3 45 + * 1 6 A 346 
+(1 6 A 356 + * i 6 A 456 * 


( 4 . 99 ) 


( 4 . 100 ) 


All of the preceding elements are sums of the constants • •*? 

weighted by A’s that are either one or zero. A possible method for computing 
© s [equation ( 4 . 91 )] by the incremental method described in paragraph 4.3 
is to initially store the 16 q's and the q Q that correspond to no failures m 
the Y registers, then when a failure occurs, change the q's in accordance with 
the preceding equations. Only the six constants j need to be 

stored to accomplish the required changes m the q's and q Q . The incremental 
computation for ar will then be similar to the uu and Aac computations de- 
scribed in paragraph 4.3* 

4.4,2 Control Compensation. Failure Monitoring, and Mode, Control 

The three compensators required for the vehicle rate loop (see Figure 
4-6) can readily be accomplished by incremental computation. The application 
of this computation technique to flight controls, for example, has been studied 
by several flight control system manufacturers, including Sperry. Due to the 
limited scope of this study, this section of the CMG control computer will not 
be investigated. 
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She failure monitoring and mode control function of the CMG control com- 
puter (Figure 4-6) generates the A 1 s for the triad conversion; the rotor 
speeds a for the steering law computer; and supplies the display signals re- 
quired to monitor the failure status of the CMG’s and sensors. The major com- 
putations required for this function are the parity check equations for the 
sensors listed in Table 4.1. A computer structure for the computations re- 
quired by the failure monitoring and mode control function is not presented 
in this study due to its limited scope. 

If. 5 CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 

The steering law computations are by far the most complex computations re- 
quired by a CMG control computer, and they outweigh the remaining computing 
functions in selecting the type of computer to employ for controlling CMG's. 
Since these computations are time-variant and nonlinear, digital computation 
is required. A general purpose computer can certainly perform these computa- 
tions but it is considerably more limited in speed than an incremental computer. 
The time required to do the steering law computations on a UNIVAC® 1819 com- 
puter (designed for airborne use) is estimated to be 6.37 milliseconds. If we 
roughly estimate that the time required to perform the total CMG control com- 
puter computations is 10 milliseconds, then a sampling rate of 100 samples per 
second saturates the computer. For some of the anticipated requirements for 
future CMG control systems, sampling rates of this order are required. It can 
therefore be concluded that the use of a central, general purpose computer, 
which is also required to perform other than CMG control functions, is not 
feasible for the type of CMG system described m this report. 

An incremental computer can be made much faster than a general purpose 
computer, but its speed depends on ihe level of serialization of its compu- 
tations. A bit-time of 0.5 microsecond is easily accomplished, and if a 300- 
bit register is used to serially store a set of 18 Y variables (such as the 
18 a’s or d’s), a computation cycle takes 0.15 millisecond. Of course, there 
are also lags in incremental computers that must be considered when comparing 
them with the general purpose type, but it is clear that incremental computers 
are much faster. Also, since computations are performed m parallel, adding 
functions to an incremental computer does not require additional computation 
time. 


In terms of complexity, it is not obvious at this stage which type of com- 
puter is the simplest. Further study would be required to determine parts 
counts. They may be quite competitive m this aspect. 
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Since incremental CMG control computers have some definite advantages over 
general purpose computers for this application, it is recommended that a further 
study be initiated to design a prototype incremental computer for a specific 
CMG control system requirement, and to compare it with a general purpose type, 
programmed for the same function, in terms of complexity (parts count), ac- 
curacy, reliability, weight, and power consumption. 
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CHAPTER 5 

SIMU DIRECTION COSINE SIMULATION PROGRAM RUNS 


5.1 DESCRIPTION OF METHOD AND RESULTS 

The purpose of the simulation of SIMU operation was to determine the error accumulation for 
large angle changes in position and also the worst-case mode of error accumulation during 
extended runs. 

5.1.1 Types of Runs 

Three types of runs were made and are described below. 

(11 Runs simulating a motion to a given angle and then back to the zero position. The 
angles used range from less than 1° to 180°. 

(21 Oscillatory runs taken near the zero reference position where the curve of data from 
the test above indicates that the worst error accumulation occurs. 

(3) Three-axis rotations where the inputs are: (al up to one radian for x , followed by 
one radian for y and then one radian for z , followed in reverse order for the 
return to zero; (b) single pulse increments in x, y and z followed successively 
until each axis has had the pulse required to read up to one radian rotation (32,768), 
then the reverse order of the same increments' of x, y and z until zero is reached. 

5.1.2 Test Results 

The results of these tests show the following: 

Figure 5-1 shows that error accumulations m direction cosine values after rotation 
to a particular angle (of 180° or less) and then returning to zero degrees were limited 
to a maximum value of approximately one part in 37. Refer also to Table 5.1. 

The error moves up steeply until approximately 2000 unidirectional input pulses have 
been accumulated. This is approximately a range of one bit of error for every two 
pulses fed the system by the time the rotatioif back to zero is completed. (See steep 
initial slope on Fig 5-1.) 

The three axis runs were made for continuous inputs of a fixed number of pulses for 
x, y and then z sequentially. They were also made for runs where the input pulses 
were increments of Ax, Ay and Az repeated in that order until the required number 
of sets of pulses were processed. Thereafter the reverse sequences of pulses were 
given to return the system to the zero reference. Table 5.2 lists the above data 
plus that for each set of two axis updates for 8, 128 and 32,768 sets of input pulses, 
the value 22.768 approximating one radian of rotation. The results show the older of 
magnitude of error for the number of pulses specified to be the same as for single 
axis operation. However, as is expected, they are not identical values. 
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FIG. 5-1 Plot of error in repeatability. 
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The oscillatory runs consisted of ten cycles from 0 to 1024 pulses and back to deter- 
mine the effect of extended testing on the maximum error rate that the system initially 
produced. CSee Fig. 5-1.) The result was that the enoi increased proportionally 
for as many input -pulses as were applied. See Table 5.3 

These tests were run to provide a reference to determine the correct operation of the SIMU. 

TABLE 5.1 

SIMU Simulation Data 



Sinale Axis 

Values After 

n(A9) Pulses 



n 


Cll 



Cl 3 


1 ,024 

037770 

000052 

165024 

000777 

165253 

027743 

2,040 

037740 

001252 

151642 

001777 

052536 

017055 

4,096 

037600 

025245 

035227 

003772 

125673 

160345 

8,192 

037002 

124477 

036241 

007725 

073552 

133023 

32,768 

021224 

050037 

135327 

032732 

124436 

043614 

106,496 

140140 

026754 

104267 

174423 

052267 

067173 


Si nal e 

Axis Return to Zero Values After n(A0) 

Pulses and 

n(-A9) 


n 


Cll 



Cl 3 

1 ,024 

040000 

000000 

037670 

000000 

000000 

041541 

2,048 

040000 

000000 

037666 

000000 

000000 

041727 

4,096 

040000 

000000 

037675 

000000 

000000 

041706 

8,192 

040000 

000000 

037675 

000000 

000000 

041662 

12,288 

040000 

000000 

037663 

000000 

000000 

041651 

16,384 

040000 

000000 

037711 

000000 

000000 

041622 

32,768 

040000 

000000 

040031 

000000 

000000 

041461 

65,536 

040000 

000000 

037631 

000000 

000000 

037303 

98,304 

040000 

000000 

037717 

000000 

000000 

037112 

106,496 

040000 

000000 

037444 

000000 

oooooo 

035560 


5.2 PROGRAM DESCRIPTION 

The program, written for the UNIVAC 1108 in Fortran IV, computes new direction cosines 
from old direction cosines using the truncation-error-free algorithm. See the flow chart 
(Fig. 5-2) to determine how this simulation was programmed. 

The nature of the iterative process is controlled by punched card inputs in the format 

described below. Provision was made for controlling iteration by selecting angular increments 

— * — * — » 

in each of the three axial senses, x, y and z . Let the three selection variables be DX, 
DY, D2 . DX=0 means skip an angular increment in this sense. DXyO means perform the appro- 
priate iteration for an angular increment in tins sense (this angular increment is assumed 
fixed). Refer to Fig. 5-3. 
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FIG. 5-2 Flow charts of simulation program. 
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computation loops for 
positive increments 


J 



\ computation loops for 
l negative increments 


J 


Choice of Angular Increment 


Choice of Angular Increment 


FIG. 5-2 Flow charts of simulation program (cont.). 
















computation for positive increments 





computation for negative increments 

'"ITER is the subroutine that actually performs the computation on the 
pair of direction cosines that are the arguments in the call statement. 

ITER itself calls another routine SCALE which scales the 3 integer 
variables of each C. to lie within 0 and 2^ - 1 ; 

The variable Z is set to +1 for positive increments and -1 for nega- 
' tive increments . 

Printing is done after each call to iter depending on the control 
- variable PR. 


Computation for New Direction Cosines 
FIG. 5-2 Flow charts of simulation program Ccont.). 
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arguments: (S,C) 

* 

S represents the three integer 
variables holding the three 16- 
bit parts of the 1st direction 
cosine, and G is used to repre- 
sent the 2nd. JPS, JEC, etc., 
represent, corresponding high order 
bits for each part. 

Computation is achieved by adding or 
subtracting the parts of S,C ,S j -2”32 i 

etc., that are of corresponding 
order. The high order bits are used to 
perform rounding and "shifting". See 
program for details. 

S'.C'.S" ,C" represent intermediary 
results which are in the same format 
as S and C . 


Actual Computation of New Direction Cosines 
FIG. 5-2 Flow charts of simulation program (cont.). 
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Scaling of 16-Bit Parts of Direction Cosine 
FIG. 5-2 Flow charts of simulation program (cont.). 


















Two further variables, L.IM1 and LB12, control the length of the iterative procedure. Two 
specific settings of LIM1 and LIM2 are of interest; 

(1) LIM1=1, LB!2 = K means iterate K times in the first selected axial sense, then K 

times in the second, etc. 

(2) LIM1 = K, LI112=1 means perform K times the sequence Cone iteration m the first 

selected sense, followed by one iteration m the second selected sense] etc. 

The program performs the "mirror-image" of the iteration sequence in order to return to 
the starting point. 

One further variable, PS , is provided to control print frequency. If PR=K , then each 
N 1 "* 1 iteration in any selected axial sense forces output of all the direction cosines. 

Input variables Cll, C12, C13 are provided to enable the user to specify different ini- 
tial direction cosines. Cll, C12, C13 are always entered and printed out as 6 digit octal 
numbers . 

The program performs only integer arithmetic, using one Fortran integer variable to hold 
the three 16-bit parts of each 48-bit Cjj , as well as intermediary results (see Table 5.4). 

Card Column 

1 

Card *1 

110 = integer, field width of 10 




Card Jf3-5 

high order 
16 bits 
c 

n 

medium order 

16 bits 

c. 

i.l 

low order 
16 bits 
c. 

i-l 


9F 

0F 

9F 


0F = octal integer, field width of F 


Cll on Card 3 
C12 on Card 4 
C13 on Card 5 



At the end of a sequence specified by one set of 
data, the program reads a new DX , DY , OZ . 

If DX‘ , = -1 the program stops, otherwise, compu- 
tation resumes with a new set of data which is 
assumed to follow card 6, etc. 


FIG. 5-3 Card input format. 
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T T r J 


r 

1 


1 


n 


^ i-Eiao'" teener 

1 "777'’ 1 77777 IKECCC 

' ~*T>T< 17777? C6"CC" 

3 "77777 177771 isOCC" 

'■) "37777 17 77 7" C^OCC" 

j "37777 17 7757 NGOC- 

"7777 1 7775? " ^CG" 1 

Y 7 7 7 7 177747 140231. 

17777 I 7 7 7 M G-lOOGG 

-i : 37777 177747 unc:'- 

-7 "’•7777 17775? r^CDDP 

-3 " ’ 7 ’ 7 7 1 77767 1 4 COLT 

-4 "7 77 ” 1 7777" 54C3C" 

-5 r ’”7’ 1 77 77 ’ 14C3CT 

-C n 27’77 *7777? C^eC" 

-7 " 3 7’ 7 7 177777 14,ia:" 

' 14033 " "Or"Q n 640CC" 



eccoo", 


DCiC"^ 

cnj"u: 

"4030“ 

uCC" 

cceeec 

J4GQ00 

GCCCC r 

no"c 

" ’777’ 

c cr r p r 

*50" oc 

o«nooc 

CCGT' 1 1 

0f10"C' 

r 77774 


c e on 03 

G4nroo 

OLCC' 1 1 

1 OO n O.'. 

"377?" 

C CCTC 0 

ooooni, 

Q4G00G 

cocr~? 

003"0" 

T 1775“ 

•JQO j-n 

cconnc 

24GGCC 

ccon*? 

uic n tr 

ri777’ 

^cc it r 

cccrrc 

ru^roc 

sesr c 7 

C“" "uC 

C 3 75 f 4 

eo »_ cj r c 

1CUCDC 

C4CCDG 

etcon^ 

iro"cc 

“176 1 t 

C C 2 3 ** p 

C0OC0O 

C4TCC0 

GGDCC4 

OCC"G" 

r» 75?" 

nonne^ 

ccocno 

G40GGC 

cccr"*' 

ncc. 

r 1761 " 

(C r, ilC 

roc cod 

04 ^roo 

ccoro 7 

cr-c"co 

r ’ 75 6 " 

— U l* / * f 1 

oCOCOo 

04 CGu 

cccr^? 

1 "C" 0 - 

77?" 

-C\“* ^ 

eocene 

G4~~3C 

CGL n ‘'2 

DCC"00 

"17744 

*- «Cou “* 

GCCGO' 


QCGC *5 1 

1 "3"G 

" 177?“ 

:c:r p 

, COOCO" 

Q4HC3C 

C CDC” 1 

D r c"o: 

“3776? 

c 

ccornt 

C4PGQG 

QGCGC‘ r ' 

1 "QH0 2 

“ 7777 - 


_ocrn 

04T'*0C 

QZGm*- 

QCC"E ", 

"" 7 77 - 


Single axis - 0 to 8 pulses and return to 0 again, 
ere 1 1 

PR»8 

SEQl 128, I-12B. 20 
S20> : 0j 2040000 
CS040000, 20. 2040800 


« 8= 

4 = 

0 = 

37 520 

tt 

37777= 

1 77740= 

40000 

= 16= 

10 = 

0 = 

352 40 

C 

37777= 

177600= 

40000 

= 2*5= 

1 4= 

0 = 

26760 

K 

37777= 

177340= 

40000 

= 32= 

20 = 

0 = 

12S00 

IX 

37777= 

1 77000= 

40000 

= 40= 

23= 

177777= 

166224 

3 

37777= 

176340= 

40003 

= 48= 

27 = 

177777= 

127754 

= 

37777= 

175600= 

400 13 

= 56= 

33= 

177777= 

55504 

xz 

37777= 

1 74740= 

40026 

= 64= 

37 = 

177776= 

I 65234 

3 

37777= 

! 74000= 

40047 

* 72= 

43= 

177776= 

54764 

E 

37777= 

1 72740= 

40100 

= 80 = 

47 = 

177775= 

122S14 

= 

37777= 

171600= 

401 43 

= 88= 

53= 

177774= 

I 44244 

= 

37777= 

1 70340= 

40222 

= 96= 

57 = 

177773= 

137774 

s 

37777= 

16 7000= 

40 320 

= 104= 

63= 

177772= 

103524 

E 

37777= 

165340= 

40 4 40 

= 112= 

67 = 

1 77771= 

1 5254 

= 

37777= 

1 63600= 

40605 

= 120= 

73= 

177767= 

73004 

3 

37777= 

161 740= 

41002 

= 128= 

77 = 

177765= 

H253? 


37777= 

160000= 

41233 

=- 8= 

73= 

17776?=' 

73014 

B 

37777= 

161 740= 

40777 

**- 16= 

67= 

177771= 

1 5274 

3 

37777= 

163600= 

40576 

*“ 24= 

63= 

177772= 

1035S4 

B 

37777= 

165340= 

40426 

=- 32= 

57 = 

1 77773= 

1 40034 

K 

37777= 

167000= 

40304 

=~ 40= 

53= 

177774= 
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K 
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0 = 
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Single 

axis - 

0 to 128 pulse 

S 

and return to 

0 aga 
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Two axis, y,z 

sequentially - 

0 to 128 pulses and return to 0 again. 
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Two axis, 

x,z 

sequentially - 0 
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to 128 pulses and return to 0 again. 
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Two axis Ax, Az incrementally - 0 to 128 pulses and return to 0 again. 
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axis 

Ax, Ay, Az 

incrementally - 

0 to 128 pulses 
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to 0 

again. 
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Two axis 

x,y sequentially 
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pulses and 
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to 0 agai n. 
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032732 

124436 

C43614 

1 04453 

101105 

0327 32 

124436 

04361 4 

1 74637 

033047 

0327 32 

174436 

04 351 4 

122024 

126362 

032732 

124436 

C43614 

1 63430 

166366 

032732 

124436 

043 61 4 

C 5 34 32 

027102 

03273? 

124436 

343614 

coocco 

G4CC27 

032732 

124436 

C436 1 4 

oaoooo 

040027 

□ 304 37 

067 HO 

046051 

000000 

040027 

025637 

174055 

1 321S1 

coococ 

04CC27 

C2256? 

035725 

C 4 33fi 

DDOOOO 

040027 

317256 

154164 

101 235 

000000 

04 OC 27 

013561 

001125 

1 7765’ 

cooooo 

G4CC27 

GC772E 

C 7 3 c 5 2 

i32srr 

000000 

040027 

003772 

125673 

1 57735 

000000 

040027 

acaooo 

OOOiQQ 

0«1 255 

32,768 pulses 

and return 

to 0 again. 


ooont 

040G0C 

000101 

0 n 0°3 p 

r. ^rn- 

oocooc 

040330 

CCCDQP 

3P11QC 

^ 4 c”r” 

oooooc 

341333 

ocoaon 

COO 10" 

” UC 3P" 

OQu pr i- 

c 4 ” ~ 

□ LOO 110 

o^c”c: 


ccocrr 

04-006 

ccccor 


r 4ccr~ 

cccroo 


cccaai 

311 , 13 ' 

r ( 1 *"''“ 

CQOODG 

04133 c 

COCCOP 

crcicc 

r4tcri 

000000 

aa^ooc 

OOOOGO 

OOG n OG 

0 4car- 

CGCDCO 

040CD6 

CCCODC 

orcicc 

T4CC r- 

C52104 

117367 

ccorcr 

crcicc 

ruccr' 

104275 

14471! 

ccgogi 

Onc”G n 

"4QC 

1 76652 

147605 

ODDCO^ 

C i ~Z'Q r ' 

~ ur cp ~ 

013613 

176300 

oocoon 

000103 

n 4o o n ” 

1 420S3 

034345 

000000 

000700 

□ 40331 

CG3722 

14*471 

OCDOOC 

GCC r CG 

r4GDrr 

013613 

176374 

00000" 

OOC n OG 

n4aor n 

1 766=2 

147743 

003001 

ooc”2: 

nuGor~ 

104725 

145375 

GCGPPO 

C r C n CC 

14 c cr- 

052! 04 

117512 

OCOCC P 

013130 

"40.“,?' 

CDODGC 

036213 

ccooor 

CCC1CC 

ruccr- 

ocoooo 

025213 

00000” 

□ G 0 '"‘G 2 

n 4aOG” 

CCCODC 

C 1 4 21 3 

ocDPor 

C"C1CC 

14 CC O' 

ooonno 

37> 21 7 

GGQCOO 

333131 

n 4CCP” 

coGmc 

C * r 2 1 3 

OLcrri 

G r cicr 

ruccr” 

000000 

C7671 5 

□GCCT” 

0 n 2 n 32 


CGOOOC 

076213 

ccdcco 

orcicc 

r40QP r 

GOOOOO 

035213 

000003 

GOO n OG 

04300- 

cocrno 

C 3 62 1 3 

ccccor 

C"C"CG 

racer” 


32,768 pulses and return to 0 again. 
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X 

1 

I 

A 'e 

32 7 6 F 

1 

IT 

C40CCC 

CCCCDC 

4C=c 

04occr 

Gcocor 

1 1 92 

040C0C 

occcor 

1 ■>293 

C4023" 

CGOOC" 

15754 

r.scooc 

OGCG09 

2 m sc 

240CCC 

003302 

2 4 5 7 5 

040000 

002000 

2 3 5 7 2 

0402 CC 

0C30GC 

32755 

C4GQC* 1 

ODQOQO 

j'»b 

0376C2 

C2E245 

-192 

r*70C2 

124477 

1 ”2 0 6 

23561= 

267 72 7 

15334 

0 34 OS ? 

050014 

*2480 

031746 

150733 

24576 

027 323 

176514 

1 S 3 3 4 

234257 

050014 

22490 

335515 

067727 

24576 

037302 

124477 

2=672 

C776Gr 

CZ 5 2 4 = 

52762 

C = OCT r 

2C-C33 

-4"°6 

- urc rr 

2C22GC 

-2192 

240CCC 

2 CC 2 2 2 

122PS 

24cccr 

2CCCC2 

163 94 

rsccrt 

CCCCCC 

2'4 3C 

240232 

CC 3302 

24576 

04 030C 

COrGGO 

25672 

Cttcoro 

ccoocc 

32766 

C4CCCC 

ccccoc 


L4rcor 

6DGOOO 

GOOODO 

csrocc 

C0C0G6 

CDDOOC 

G4CC0C 

GDCOOO 

□00000 

046*00** 

OOCODO 

000000 

04000** 

OOOCCO 

OOOGCG 

6.4020*7 

G 000 3 0 

000000 

G403C0 

OOCOOD 

oooooc 

0403 00 

000300 

cooooc 

C4Q300 

000030 

00DC0C 

63=227 

CGCC6C 

CCGCCC 

C 36 24 1 

CDCCCC 

coocoo 

C 7 2 35 7 

C0CC6.C 

cooocc 

175554 

200000 

oocooo 

0S4306 

GOOOCO 

000000 

172114 

000000 

□ 00000 

176151 

oocooo 

oooooo 

373033 

□03300 

0000 DO 

G36G23 

000000 

oooocn 

635552 

COGGCO 

60D00G 

C4016 3 

CDCCCC 

cooccc 

046163 

COCGGP 

GCDGCO 

C4C163 

CDCCCC 

cooccc 

G4 C lb 7 

COG OCO 

coccoo 

CS0163 

COCCGP 

cooooo 

C42153 

00C3CC 

ooooco 

C4C163 

0OG03C 

COOCG3 

0431 53 

GOCOCQ 

GDDOOC 

640163 

GOGCCQ 

coocoo 


C4GCCC 

CLCCCC 

060=00 

C4CTGC 

CCCCOC 

C n G”CD 

04rccc 

CCCOCf' 

crcoo 

340330 

□ Conor* 

□"□"DO 

043000 

230C2 r ' 

ODO n GG 

0400 Ob 

61C0Q0** 

OOCOOO 

CU^OOC 

2G0C03 

0n2 n G0 

C4030Q 

GCOC *** 

Dro-oc 

040200 

OuOCOO 

QQQ-QO 

G4rccc 

CC3772 

1 ? 5= 7 3 

c4rcoc 

CC7725 

673=52 

CC-CDG 

C 1 3561 

CC 1 1 25 

040030 

01 7256 

164164 

040020 

02256? 

□35725 

04n003 

025637 

1 74^55 

042027 

01 7256 

154164 

043327 

013561 

001 1 25 

042C27 

007725 

073552 

G40C27 

CC3772 

125= 7J 

C4CC27 

cccccr 

GCC-CG 

C46C2 7 

Gtocor 

CfC"CC 

C4CC27 

COCCC" 

COC"CC 

C4CG27 

CCCGOC 

crc'cc 

C4CC27 

CCDCCC 

cgctc 

0 4 32 27 

0C3CC*- 

orc-oo 

040327 

0CC022 

□oc ri ac 

J40D27 

OOOCCr 

ogc-gc 

C4C027 

CCOODO 

GCC r C 0 


rucof 

64cc-- 

64cor , 

C4C20 

me or' 

Gaaor- 

rooor* 

0400-' 

04020" 

1=634= 

1 ** 0 ?’ 
1*61?* 
1 n l 44' 
G 43 274 
1 7225* 
I'M 22 = 
1 * 7 65 * 
1 323C" 
1 = 77 3 = 
C 4 12 6 = 
C 4 1 2 5 7 
64)28** 

C4 12 5 c 
C 4 12 5 = 
“4! 25' 
r 4i 25= 
r 41 25 1 
64125= 


Three axis 


x,y,z sequentially 


0 to 32,768 pulses and return to 0 again. 


1 1 r 


4095 1 3276P 


G 

040000 

GOOOOO 

□40000 

OOOOOO 

OOOOOO 

04 GOOD 

ooocor* 

000=00 

0 40 00 = 

4G9S 

0 3 7600 

063451 

□0601G 

000177 

121314 

047357 

OG 3 764 

155427 

1 73 57 = 

*096 

0376CI0 

05?476 

153331 

000177 

1 21314 

047357 

003765 

0SSP30 

1 30 74* 

3 19 2 

037005 

06737S 

GG324C 

000772 

1 20254 

1 50745 

00765? 

134*73 

1 **2 33 5 

SI 92 

0 37CG5 

047556 

0744 37 

000772 

1 20254 

150745 

OC 765 3 

0*2*GS 

041 sr* 

172 33 

035632 

166516 

0172 77 

002145 

024625 

17*041 

01 334* 

111 422 

C*1 067 

IT'S 

035632 

1 376CG 

1G=*S45 

CD2145 

G246Z5 

1 7 ? C4 1 

013344 

CC51G7 

1474 r7 

15 334 

0341 24 

G23505 

070720 

003653 

1 73037 

153720 

Cl 6S45 

144*03 

061 4*7 

15334 

0341 23 

166172 

06664? 

003653 

1 73037 

153720 

016546 

035153 

0 73 3? 7 

26430 

C3 2 1 1"? 

G4C642 

147771 

605664 

160533 

Cl 662 6 

021 375 

D 66 ?C 6 

17277= 

2C43Q 

032112 

175647 

105725 

005 66 4 

160533 

015626 

021375 

1531 35 

0 31 51 * 

24576 

027646 

06 7021 

044715 

010137 

1 34555 

0127 32 

023575 

157*10 

C3312S 

24 = 76 

927640 

017425 

C2537S 

010137 

1 34555 

012732 

G23576 

07 7 4 1C 

1 41 57* 

2-672 

0251 7** 

13512* 

12?4 35 

012607 

□40163 

007237 

025303 

032-53 

n ll 77* 

25572 

CP517Q 

112311 

153765 

0 1260 7 

040163 

007237 

025303 

124444 

1 11 425 

32758 

02237= 

153410 

150573 

01S4O2 

052652 

177760 

, 026263 

060*21 

052171 

3:?6e 

622375 

67664 1 

15173= 

C1E4D2 

052652 

17776C 

C2S26* 

115714 

17444* 

-4,-55 

67517*- 

11231 1 

153732 

012 = 37 

112771 

172752 

025303 

06 7°25 

1 3626' 

-4696 

6?= 1 7= 

1 1231 1 

16373? 

Cl 2607 

C40163 

C*'71SC 

025367 

l r 4 4 4 4 

1114*7 

-3192 

027640 

Cl 7425 

G25342 

010140 

0041S1 

07 1 275 

023576 

Cl 71 1 1 

**20 605 

-3192 

02764C 

017425 

025342 

u 1 0 1 3 7 

1 34555 

012554 

023576 

037410 

1 41 534 

-12288 

622112 

175647 

105732 

CO 5 6 65 

623526 

1 ? 1 ? 1 3 

021*75 

1*7*63 

6*51 r? 

-12288 

-3211? 

1 7SS47 

1067 32 

005554 

160533 

015473 

021375 

153135 

** 31 36* 

-15384 

2341 23 

I6G17* 

u 565C7 

003654 

0 30354 

0*5 356 

01 6546 

0*5423 

"4753' 

- 153S4 

034123 

16517? 

0566 07 

003653 

1 73037 

153560 

016546 

035153 

"73244 

-22430 

035632 

137606 

130316 

0 021*4 5 

053S3S 

177642 

01 3344 

090 = 7 S 

247045 

- 22 4 8 C 

035632 

1 3 7606 

100516 

002145 

024625 

171736 

01 3344 

005107 

1 47 43! 

-24576 

037005 

047650 

G74463 

000772 

1 40003 

033424 

□07653 

070 *20 

1 61 345 

-24576 

C37C05 

C 4 7 650 

674463 

000772 

320254 

1 5 C 5 7 5 

GG7653 

0 * 2 *D 6 

641625 

-23672 

03750C 

C5->476 

153311 

000177 

1 31267 

Q0 n 6 6 5 

003755 

0544 31 

056102 

-29672 

6*7606 

CS2476 

15*313 

CGC177 

121314 

C472G3 

GC 3 76 = 

C=5-*C 

1*67== 

- 32763 

040900 

OOOOOC 

C373C7 

GOOOCC 

OOOOOO 

037675 

OCOCOO 

COOCOO 

- 42 71 7 

-32763 

040C00 

00000 " 

G 3 7 6 07 

OOOOOO 

oooooo 

037S7S 

OCOCOC 

000 n Q 0 

C«2 71 7 

Two 

axis 

Ax, Ay incrementally - 0 to 

32,768 pulses 

and return 

to 0 again. 


TABLE 5.2 (cont.) 


108 



C J ' 1 


4C96 

□ 

1 

040000 

32758 

2G0C0C 

040000 

oooooc 

000000 

340000 

000009 

0Q0-3G 

C400 r ~ 

4096 

0374cr 

135147 

121541 

1 7401 3 

022350 

104712 

003765 

055" 30 

1 3C214 

4396 

237400 

125175 

057460 

174012 

1 23347 

322345 

GO 3765 

055"3C 

1 30214 

3192 

03601? 

137246 

120273 

170125 

043504 

176160 

007653 

032306 

041163 

3192 

C 36 01 2 

117520 

133456 

170124 

1 47457 

117171 

007653 

032306 

C41 1 63 

12293 

C3 3465 

1 26 3 2 4 

9505 54 

164434 

063356 

057443 

01 3344 

D95137 

1 46735 

12288 

233465 

277414 

146345 

164433 

1 77203 

031503 

31 3344 

CC5 1 0 7 

1 46735 

15339 

03025'' 

G33 7 C n 

120444 

161232 

C 33075 

0171 35 

016545 

C351 53 

C 72 67 3 

16384 

C 30 247 

154364 

126115 

161231 

1 52355 

03C6SG 

016546 

035153 

0 72 67 7 

2043G 

024225 

G3S512 

016777 

1 56402 

1 11071 

105501 

021375 

153135 

03113" 

20 4 90 

024225 

I73S1S 

170506 

1 584 0 2 

040415 

053250 

0 2 1 37 S 

1531 35 

031 13" 

24576 

017500 

1 OS 4 4 6 

C 334 14 

1 S42C2 

020070 

045234 

023576 

C 37 4 1 0 

1 43 345 

24576 

017500 

03705 7 

234373 

1 54221 

1 60667 

077514 

023575 

0 37 4 1 C 

1 41 34 = 

23572 

012361 

Q7743 7 

0474 42 

152474 

1 457 IS 

366255 

025303 

1 n 4 4 4 4 

1 1 1 361 

23572 

012361 

024623 

136332 

152474 

1 207 52 

141777 

025303 

104444 

1 1 1 361 

32768 

004777 

C52 2.52 

073662 

IflS 14 

12^X57 
1 1 50 70 

026042 

C26263 

115 7 14 

1244 71 

32768 

004772 

175503 

140032 

151514 

156054 

026263 

115714 

1 ?4 471 

-4096 

Cl 2360 

152014 

100323 

152474 

1 20752 

145044 

C253C3 

131406 

110 3 r 4 

-4 0 96 

Cl 2 3 6 1 

224 6 2 2 

1 36376 

152474 

120752 

142044 

G 2 5 3 C 3 

1C4U44 

11144! 

-SI 92 

0 1 7477 

167455 

075757 

1 54201 

1 60667 

D77637 

023576 

076510 

1 70234 

-3132 

01 75CC 

G37CS? 

0 34 2 GO 

1 54201 

160667 

077637 

023576 

037410 

1 41 507 

-12286 

024225 

13CS22 

171332 

156402 

C60415 

C 5 7 406 

021376 

C 2 7C 10 

1 = 55 7 7 

-12288 

024225 

173516 

170277 

1 56402 

040415 

053406 

021375 

1531 3S 

031 273 

-15384 

C302«7 

117047 

1527 36 

161231 

1 523 55 

033736 

016546 

J 1 5= 72 

1 6653" 

-16384 

C 202 4 7 

154364 

126CC3 

161231 

152355 

C3C736 

C lb 5 96 

C 351 S3 

073074 

-20430 

033465 

G505Q4 

044723 

164433 

1 77203 

031575 

013344 

074 ■’62 

1 1 7405 

-2048C 

033465 

077414 

146316 

1 64433 

1 77203 

031575 

013344 

005 137 

14726= 

-24576 

C 7 6 C 12 

C7777 1 

1 52725 

170124 

147457 

117242 

GD7652 

1 2£ 7 3 2 

06123 1 

-24576 - 

C 350 1 2 

117520 

133566 

1 70124 

1 47457 

117242 

00765 1 

0323Q& 

041 47 7 

-23672 

277400 

115222 

2365 32 

1 7401 2 

1 23347 

022572 

00 3 7 6 5 

154 n 31 

17323" 

-58672 

C77 4 or 

125 175 

067613 

174CL2 

123347 

C22572 

CC37SF 

0 55" 3C 

130611 

-32758 

037777 

177777 

137772 

OODOOO 

OOOOOC 

036571 

OOOOOC 

103"00 

041 77<- 

-32768 

C4000C 

QCOCOO 

037772 

000000 

OOOOOC 

036571 

000000 

000"00 

C 4 1 77 r 


Two axis Ay,Az incrementally - 0 to 32,768 pulses and return to 0 again. 


i n i 

3 a 1 '""“l 2 


c 

u r f r 

cor: o'" 

L = " J_r 

GOCLO 

ccc r oc 

'frr' 

oco"or 

r r ^ r 

m 11 " 


r 7 75rrr 

lo?451 

: 7 c 

1 7471 7 

C 221 

1"4 1 1 7 

Q001 77 

1 ?1 *1 4 

- 4?3«K 


" 7 75G r 

95?4 7 5 

153 : 1 = 

1* 7 4 3 1 7 

1 22747 

14^747 

390177 

i?i 7 i 4 

" 4 7 354 

1 12 

> 

06737= 

053215 

1 7012= 

0 43504 

17^254 

CCG77 7 

, 1 7 " 7 5 4 

1 5C 76' 

! 3 ■ 

r 77 -rr 

947=5" 

c 7441 = 

17C 1 24 

14 = 47^ 

C 7 = 106 

C9G77? 

1 "9 7 5 4 

7 = F 76 = 

1 ? =3 


IS*! 1 0 

1 7 2 4 7 

1 b44 34 

G 

2 " ~ ^ 3 ? 

092145 

3^4'35 

17»;c/ 

•J 
1 J 

(7 

~ 7 5 5 3 ■> 

127335 

i:"5i5 

1 £4433 

1 7? = 7o 

119311 

03214= 

3*>«i 

1 7 2 7 = 4 

1 j Z 5 4 



* 1 ~ j 7 

1 61 23 7 

033075 

2 ' - ZZ Z 

3 9 365 7 

1 73" 3 ? 

I = 3 74 " 

iz :s4 

" 7 412 7 

1EC17" 

c c f : c c 

161 231 

14 26 25 

:rtni 3 2 

9 C 3 6 5 7 

1 7 3" 3 7 

1=274" 


0 3 ? 11 7 

G4CS'4? 

1477 

1 5640? 

1 1 1071 

n4"4i 


150=33 

“ 1 65^^ 

Z^UPC 

0321 1 ? 

17=547 

1 " = £6 7 

156402 

024643 

046027 

005664 

160= 77 

01656" 


0 "75 4? 

057321 

rmijv" 

1 S42C 2 

0 209 70 

394525 

Cl 01 7 7 

1 74=55 

r 1 ? 77 7 

24 = 73 

027649 

01742 = 

02=330 

I 54201 

1 40367 

1 7 6 I 5 7 

01 31 3 7 

174=55 

" 10 77’ 

-4: 9o 

n*si 7- 

1123 1 1 

1 5 7 5 ” 1 

1 52474 

C 73333 

1=6333 

01 Z5C” 

049 163 

7 t ? -, 

- -i°z 

~ - 7 s 4 " 

317425 

C" = 71 

15 4*201 

160 F 6 7 

3 7 7 15 1 

C1C14" 

0 r 4 T 5 1 

rail’s 

- 1 1 ^2 

r 7 7640 

037425 

.-’=>71 

1 54201 

1 40367 

l 7=223 

010137 

174=55 

0 12 £l« 

-1 2298 

03 7 1 ! 7 

17S547 

105553 

1 56402 

04041a 

352537 

GO 566 5 

0?3=2a 

1 21 24 = 


I 1 ? 

1756- 7 

1 C?S 6 ^ 

1564C2 

C24S43 

1 4 6 3 5 b 

C. 55664 

160=23 

r 1 S 5 ?f 

- 1 5 7 7 4 

“341 7 " 

1351 7" 

' 4 7 

16123 1 

1 52355' 

JVMC2 

303654 

0 ’G 7 54 

~ 354 r- 

-16144 

"!4l^ 

165 1 7 7 

w50!) 4 3 

« bl 231 

1 42=?a 

2 ''ft 4 60 

0 r 3 6 5 1 

1 73"37 

1 5361 - 

- t =* * r 

r 7 s= 

13750= 

100454 

164433 

17 720 3 

G 7 "65C 

CL? 1 4 = 

0=3=35 

1 776F = 

- ?Oh 97 


137505 

lb"Hi4 

1 5443 3 

1 7267„ 

1 ^773 

□02145 

074=25 

1 71 751 

- 24 = 7:5 

~ 7 7 f '2 c 

04735" 

0744 30 

1 7C124 

1 4 7457 

11=352 

00077" 

1 4u"C 2 

" 3344 ’ 

“?4 C 7 b 

7 r r- 

C 476 5 r 

C 7 4 4 3 ^ 

170124 

145472 

C’f"7£ 

C CO 7 7" 

1 "9"5 4 

1 =CS 1 4 

- 73 £7 ? 


05 _, 4 7 5 

i 5 ^ fc. 39 

1 74012 

1 233 = 7 

""1512 

3C9 1 7 7 

1 71 "57 

57" 

- 7 * :» 72 


C5-^4 7 ■; 

! 3 7 ca 6 

1 74012 

1 22747 

1 4"33£ 

3901 77 

1 7] 7 1 4 

"4721 ! 

"Ji 7f» fc 


Z 2C1Z~ 

C 7 7 27 7 

CCCC"" 

coccor 

C 7 " 790 

occcor 

C 0 L"CC 

" 7 766 = 

-: 77 58 

'’40000 

oocao" 

o’ 75 77 

Gera-" 

OODOOO 

374 793 

occrc" 

30u"o3 

" 7 7 £ 6 ' 


Txvo axis Ax.Az incrementally - 0 to 32,768 pulses and return to 0 again. 
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1 

l 

i 




i 

2276’ 



f 

rurorr 

CCCC2C 

C4CCC0 

COOOCP 


r?74r l 

017531 

14323? 

179217 

<,196 

0 274 Cl 

007173 

037401 

1 7421 7 

4C96 

0 3 7 4 Q- 

177632 

. 115160 

1 7421 7 

9192 

076020 

022717 

17 7 7S« 

171167 

= 132 

036020 

0CC531 

C04273 

171187 

3192 

03601 7 

163237 

014117 

1 71166 

l 22 7 5 

033520 

034055 

031073 

1 66777 

l’ZE3 

C33S2C 

CC1 3 1 4 

1 11036 

166777 

I 22EE 

C33S17 

15731? 

166147 

166776 

1S3E4 

C7C772 

C2C225 

C 35 4 1 3 

165535 

J *2 B4 

030371 

1545C3 

175326 

16SS35 

163E4 

0207 7 1 

1 27777 

00403" 

165535 

2?4FC 

074531 

115442 

147531 

165261 

22 it 3C 

- C24531 

041470 

063422 

1 6526 I 

22930 

024431 

C 1 4 2 31 

074112 

165261 

24576 

020311 

151363 

024561 

1 5600 3 

29 573 

020311 

G66214 

074165 

166003 

29576 

020311 

042222 

137142 

1 6600 3 

21S72 

013561 

01 7930 

0237 45 

1 67S02 

2E672 

[■13660 

126377 

037020 

167502 

2ES72 

C1366r 

1C56C4 

163C73 

167502 

22766 

0071 7t 

035512 

C722CC 

172111 

3:788 

00717’ 

14C726 

011367 

172111 

5? 78 £ 

007 1 7 1 

l'24 5 = 1 

027701 

172111 

-4696 

01 2'f r 

C 1’4 55 2 

050222 

16 ? 5 03 

-4036 

013=6* 

105604 

163035 

16 7 5C? 

-6C=>6 

013=60 

1056C4 

163036 

167502 

~?1°2 

020714 

1 5705 3 

065151 

166003 

-31 32 

020311 

042222 

157176 

1 6600 3 

-31 32 

C20 311 

242222 

1371 76 

1 66003 

-12293 

n? 4 = =- 

140260 

0125 21 

155261 

-12283 

224531 

014233 

074106 

1 65261 

-12233 

024531 

014233 

0741 OS 

165261 

-15394 

030371 

2642S4 

13055? 

1 65535 

-IF 364 

072’ 7 3 

1Z7777 

004 121 

165535 

- 152E4 

r 7 o 7 7i 

127777 

004 121 

165535 

-:046C 

023517 

12455C 

C23E9 I 

166777 

- 2 rusc 

033617 

157312 

166402 

166777 

-204SC 

C77517 

157312 

166402 

166776 

-24576 

C’&m 

141624 

OC 266 1 

171166 

-24576 

"3S"17 

1 6 3207 

014516 

171166 

-24576 

"’SCI 7 

163207 

014516 

171166 

-23572 

027402 

167277 

17=000 

1 7421 7 

-23672 

037400 

177632 

115590. 

174217 

-23672 

037400 

177632 

115540 

1 74217 

-32768 

037777 

177777 

140211 

OODQOO 

-22 7 66 

24crcr 

ccccco 

0 4 C 2 1 1 

CCCOCC 

-2276S 

r4ccc" 

OCOGCC 

C 4 0 2 1 1 

CCCOCC 

Three 

axis 

Ax , Ay , Az 

incrementally - 


coocoo 

0 9 = t . G C 

COOCOO 

crc"cc 

r4top= 

16 6 27 C 

157502 

0C41 S’ 

CO 2 ’ 1 4 

ri 06 r= 

1 66270 

157 502 

004157 

% 101316 

C 3 7 53’ 

067266 

1=0455 

0C«1 57 

101 ’16 

0 7753’ 

041732 

u 4 8 0 4 2 

010570 

1 35’14 

1 472SS 

C41732 

045042 

010571 

G’l ’54 

- I 72 3?" 

145672 

06 200 2 

010571 

031 ’54 

1 7237" 

062164 

12=334 

01 5260 

114776 

1 =676= 

062 164 

12=334 

Cl 5261 

GO 2 = 37 

0143=7 

172774 

144505 

015261 

003537 

014357 

133633 

C2ES55 

02165C 

C67“4C 

16756" 

133633 

C25S55 

C2 1 6 5r 

190=24 

1 F 46 1" 

052647 

171C53 

021650 

150=24 

16451’ 

130765 

06354 7 

02S765 

CC5’22 

042574 

1 30765 

053547 

025755 

055=05 

C 21 72’ 

057503 

005524 

025765 

056=05 

021 727 

055677 

067767 

031463 

033=64 

027061 

G 556 77 

057757 

031 463 

074=07 

066661 

C 1 5054 

137331 

031453 

074=07 

C 66 661 

1 70306 

172684 

039414 

051=71 

1 71 367 

IJ0306 

172664 

C 3 4 4 1 4 

110=33 

127606 

140545 

136460 

034414 

110632 

13 76GF 

116341 x 

C36764 

C3&97 7 

171UC 

1=3037 

116341 

02 F 764 

03647’ 

l ’ 7 = CO 

151G77 

077751 

151665 

03647" 

137=00 

1 5 1G 7 7 

031=77 

CCCE52 

034419 

1=1166 

06 75 F 4 

031577 

000652 

034414 

131425 

1 4 52 Ou 

140545 

13E371 

034414 

110=33 

13755" 

ICO 2 2 3 

154447 

031463 

1F1’22 

1726 00 

1 00223 

1S9447 

031463 

132=00 

l 51 3r" 

015054 ’ 

137237 

021482 

074=07 

"5660’ 

1 33455 

170440 

02575= 

1 = 5 ’2 3 

1 8327= 

1 33455 

1*70440 

025755 

1 2 4"4 2 

036511 

0 57503 

005571 

02576= 

056 =05 

021 671 

116371 

066 1 21 

021651 

056115 

007558 

1 1 6 ’ 7 1 

066121 

021650 

175131 

173307 

C52647 

120343 

021650 

1 =0924 

16454? 

C25466 

175243 

D15261 

1 1 5"CC 

IPGCC1 

025466 

17=243 

D 1 526 1 

G 25 = 4C 

17371= 

172724 

143744 

C 1 526 1 

0=3=37 

C 1 43 3 1 

167254 

161543 

C1D571 

1 4 2 7 36 

152134 

1 67254 

161543 

013571 

046=77 

C ’5 1 

1 45672 

061 357 

010571 

G 3 1 7 54 

1 7225" 

0 77675 

162334 

304150 

007561 

0 37 36 = 

077625 

162334 

004157 

110=57 

0 5024" 

067266 

150143 

004157 

1C1 7 16 

C 3 7 352 

OODOOD 

037025 

0C0C30 

100"00 

037137 

COOOOC 

037025 

CCGGCO 

crc’cc 

037137 

COOCOO 

C3 7 025 

cccccr 

crc n cc 

C 3 71 3 7 


0 to 32,768 pulses and return to 0 again. 


TABLE 5.2 (cont.) 



1 

1 

I 



1 

1 

'9 

1 


□ 

040000 

ooooon 

040000 

ooooon 

1' 

037777 

177777 

140000 

000000 

H 

037777 

177777 

0400 DO 

I 77777 

2, 

037777 

177777 

040000 

• 177777 

2, 

077777 

177774 

137776 

177777 

2. 

037777 

177774 

040000 

1 77777 

3 

037777 

177774 

040000 

177777 

3 

0 377177 

177771 

13777? 

177777 

3 

0377177 

177767 

040000 

177776 

4 

037777 

177767 

040000 

1 777:76 

4 

037777 

177763 

137764 

177776 

4 

0 7 77 7-7 

17 7760 

040000 

1 77776 

5 

037777 

177760 

040000 

177776 

5 

037777 

177753 

137754 

177776 

5 

037777 

177747 

04000n 

1 77775 

6 

0377777 

177747 

040000 

1 77775 

6 

037777 

1777:41 

13774? 

1 77775 

6 

077777 

177734 

D400DC 

17T177S 

7 

0 3 7 7 7*7 

177734 

040000 

1 77775 

7 

037777 

177735 

137726 

1 7 77:75 

7 

037777 

177717 

OOPOOO 

177(774 

a 

037777 

177717 

040000 

177774 

a 

037777 

177707 

137710 

177774 

a 

037777 

17 7700 

C40D00 

1 77774 

9 

037777 

17 77-00 

040000 

1 77774 

-l 

7737 7 7 r 

177700 

040000 

1 77774 

-2 - 

037777 

177707 

137710 

1 777:74 

-2 

037777 

17 77*1 7 

OMOoon 

1 77774 

-2 

037777 

177717 

040000 

1 77774 

-3 

0 3 7 717 7 

1 77 72 5 

1 37726 

177775 

-3 

037777 

177734 

040000 

177775 

-3 

0377(77 

177734 

040000 

177,775 

-4 

037777 

177741 

1 37742 

177775 

-4 

037777 

177747 

040000 

1 77775 

-4 

037777 

177747 

C4C00C 

177775 

-s 

037777 

177753 

137754 

1 77776 

-s 

037 7(77 

1 77 760 

0 4000 0 

177776 

-5 

037777 

17 7760 

040000 

1 77776 

-6 

037777 

17 77 S3 

1 37764 

1 77775 

-6 

037777 

1771767 

040000 

177776 

- 6 

037777 

177787 

04nooc 

1 77776 

-7 

037777 

177771 

13777? 

177777 

-7 

037777 

177774 

040000 

1 77777 

-7 

03777*7 

177774 

040000 

1 77777 

-8 

0377*77 

177775 

137776 

177777 

-8 

03777-7 

17 7777 

040000 

1 77777 

-8 

037777 

177T77 

040000 

177777 

-9 

C 3 777 7 

177777 

140000 

000000 

-3 

040000 

00 oooo 

040000 

000000 


ooooon 

040000 

DOOOOO 

oooooo 

040004 

oooooo 

0400DD 

ooooon 

loonoo 

037777 

1 00000 

040001 

000000 

1001QO 

037777 

1 00001 

040001 

000000 

100001 

n 3 7 774 

1 00001 

040001 

000001 

OOOOOI 

037771 

000001 

040007 

oooom 

OOOOOI 

n 3 7 77 1 

000003 

040013 

000001 

000003 

0 37 76 5 

000003 

040013 

oeoooi 

10D n 03 

03775? 

100003 

040032 

0000 E>I 

100003 

03775? 

100006 

□ 401042 

00 0001 

1 00*106 

037740 

100006 

040042 

C00002 

C0OOO6 

077712 

000006 

040076 

000002 

000006 

037 71 ? 

000012 

040116 

000002 

00Q n 12 

03767? 

000012 

040116 

0000072 

ioori 2 

03 76 25 

100012 

040173 

000007 

100012 

037625 

100017 

040223 

000002 

100017 

037573 

100017 

040223 

000003 

DQ0 n 1 7 

□17503 

COO 01 7 

040325 

0000073 

oc or 17 

r3?sn? 

000025 

0401371 

000003 

000^25 

03743? 

000025 

040371 

000003 

100025 

037320 

100025 

040524 

000003 

100025 

037320 

100034 

040604 

000003 

100034 

037236 

100034 

040604 

000004 

000034 

037064 

D0D034 

'040774 

ocooou 

003034 

037064 

000044 

041074 

000004 

O n Q044 

036764 

000034 . 

040774 

000004 

000034 

037064 

1 00034 

040603 

000004 

O0C n 34 

037064 

1 00034 

040603 

00000"' 

100034 

037235 

100025 

040521 

000003 

10O"25 

03731 5 

C00025 

04G36S 

ococc 

ICCC25 

037315 

000025 

040365 

000003 

000025 

037433 

C00017 

040321 

000003 

000017 

037477 

100017 

040216 

OCODD0 

CCC017 

T374 77 

10Q017 

040216 

000002 

10D n l 7 

037565 

100012 

040164 

CCCCD? 

1 CC 1 * 1 Z 

037616 

000012 

040105 

000002 

10001 2 

0 3761 6 

OO0O1 2 

040106 

000000 

OCQ0 1 2 

03766? 

000006 

040066 

000002 

000 n 06 

03770? 

100006 

040031 

000002 

OPQ0Q6 

03770? 

100006 

0400 3 1 

OCOOOl 

1CC0Q6 

C 377 27 

100003 

04001 7 

OOOOOI 

lG0 n Q3 

03773? 

000003 

037777 

OODOOl 

inc'’D3 

T377 37 

000003 

037777 

OOOOOI 

000003 

037751 

onoooi 

037773 

000001 

000001 

0 3 7 755 

100001 

037764 

OOOOOI 

QOG0O1 

037755 

100001 

037764 

ooooon 

100001 

037760 

100000 

037762 

OOOOOO 

inoooo 

037760 

oooooo 

037763 

OOOOOO 

100000 

037760 

oooooo 

037760 

oooooo 

000000 

037760 


Three axis Ay,Az t Ax incrementally - 0 to 8 pulses and return to 0 again, 

TABLE 5,2 (cont.) 


Ill 



U> U5 


i 

1 

0 

1 

2 
2 
Z 
3 
3 

3 

4 
4 

4 

5 
5 

5 
£ 
£ 

6 
7 
7 
7 
a 
e 
a 


-1 
-1 
-2 
-2 
-2 
-3 
-3 
-3 
-it 
- 4 
-<( 
-s 
-5 
-5 
-6 


-T 

-7 

-7 

-8 

-8 

-8 

-3 


1 I 


1 

- 9 

1 

D4OOC0 

COOOOO 

C40000 

037777 

ff7777 

lWOQCD 

037777 

177777 

140000 

037777 

177777 

037776 

037777 

177775 

137776 

077777 

177775 

137776 

077777 

177774 

037770 

037777 

1 7717 71 

137772 

037777 

177771 

1 3777? 

037777 

1 77787 

037756 

037777 

177783 

1 377 64 

037777 

177783 

137784 

037777 

17 7760 

037740 

□37777 

177753 

137754 

037777 

177753 

137754 

037777 

17 77 47 

037716 

037777 

177741 

137742 

037777 

177741 

137742 

037777 

177734 

037670 

0377177 

177725 

137728 

037777 

177725 

1377 28 

C37777 

177717 

037636 

037777 

17 77X7 

137710 

037777 

‘ 177707 

137710 

C 3 77 7 7 

177700 

□'3 7 60 0 

03777 7 

17 7T07 

137710 

0 3 7 7>7 7 

1 77707 

1 3771 D 

037777 

17 77-17 

0376 36 

037777 

177725 

137726 

0377.77 

177725 

137726 

037777 

177734 

037670 

037 7t7 7 

177791 

137742 

07 77 77 

177741 

137742 

0377-7 7 

177747 

037716 

077777 

177753 

1377 54 

C377i77 

177753 

137754 

C377TT 

177760 

037740 

037777 

177763 

137764 

C37777 

177763 

137764 

037777 

177767 

037756 

0377J7 7 

177771 

131772 

037777 

177771 

13777? 

0377177 

177774 

037770 

037777 

177775 

137776 

037777 

177775 

137776 

037777 

177777 

037776 

037777 

17 7777 

140000 

037777 

177777 

140000 

OOODCO 

OOOODO 

040000 


JiQOOOO 

177777 

JJPOOOD 

,040000 

JOOOOC 

□40000 

177t777 

100000 

040000 

177777 

1 00000 

040000 

177.7 77 

000 OOD 

C40C03 

177.7 77 

000001 

04C0D7 

1 77777 

000001 

040007 

1771776 

X00001 

040021 

1 77776 

1 00003 

040031 

1771776 

100003 

040031 

1 77776 

000003 

040056 

1 777176 

000005 

040076 

1 77776 

000006 

040076 

177775 

l ooooe 

040142 

177775 

100D12 

040172 

177775 

1 00012 

040172 

1771775 

000012 

0«0261 

17 717 75 

DOOO 1 7 

040325 

1 77775 

000017 

040325 

17 717 74 

1D0P17 

040443 

1 77774 

100025 

000523 

1771774 

100025 

000523 

1 T7774 

000025 

000674 

1 77774 

000034 

^ 000774 

177>774 - 

000034 

040774 

177774 

000034 

040774 

177774 

000025 

000674 

1 77774 

100025 

040522 

177774 

100025 

040522 

177774 

100C17 

CB0440 

177775 

000017 

000321 

1771775 

000017 

Q0O321 

1 77775 

000012 

040255 

177775 

100012 

04C165 

1 77775 

100012 

040165 

1771775 

100006 

000133 

177T76 

000006 

040066 

1771776 

000006 

040066 

17 71776 

000003 

□40046 

177776 

100003 

040020 

17 71776 

100003 

040020 

1777:76 

100001 

040006 

17 717 77 

000001 

037773 

177777 

000001 

037773 

177777 

000000 

037767 

1771777 

IDOGCO 

037763 

1 77777 

100000 

037763 

177777 

100000 

037761 

COOOOO 

COOOOO 

037750 


000000 

GGOIOO 

04000* 

OOOODO 

ocoroo 

r4*D0 0C 

DDOCOO 

oro*oi 

040000 

000001 

400*01 

037776 

000000 

irDnci 

037776 

QDDOOO 

100*03 

r,377 74 

□00001 

000*03 

037765 

OOCOOl 

oro*Q3 

037764 

000001 

000*06 

□ 37 761 

G00CC1 

100106 

037743 

OOCOOl" 

inonos 

nl774 1 

OCOOOi 

100*12 

0 3 7 72 1 

000002 

□00*12 

03767? 

000002 

010*1 2 

037673 

ocooo? 

000*17 

T37SS? 

000002 

100*1 7 

137574 

000002 

100017 

037574 

oooco? 

100025 

03754? 

OCOOO 3 

000*25 

037437 

000003 

000*25 

03743? 

000003 

000*34 

□37373 

0G0003 

100r34 

C37237 

0C0003 

100*34 

03723? 

QCGC03 

100*44 

037155 

0C0009 

G00044 

036764 

0CDC03 

100*44 

03715# 

000003 

1 COO 34 

037234 

000003 

100*34 

037234 

CCDC03 

000034 

03736? 

000003 

DCCP2S 

0374 37 

000003 

000*25 

137433 

OGOOC? 

1C0025 

0375 38 

000007 

100*17 

037568 

OCOOO? 

100*17 

T3 75 65 

000002 

000*17 

13764? 

OGOGO? 

0CGri2 

C3766? 

000002 

000*12 

03766? 

QDOOOI 

I CG* 1 2 

037716 

ooooo-i 

100*06 

037726 

0000011 

10G1Q6 

037 726 

OCOOO 1 

000006 

0377 45 

0C0C31 

000*03 

037751 

OOCOCVl 

GCD*Q3 

037751 

OOOODO 

100*03 

037757 

oooooo 

100*01 

*77757 

COOOOO 

1CG*01 

037757 

ODOOOO 

000*01 

037761 

OOOOOO 

000*00 

03776* 

OOOOOO 

□OCOGO 

037760 


Three axis Az,Ax,Ay incrementally 


0 to 8 


pulses and return to 0 again. 


TABLE 5.2 (cont.) 


112 



I 1 1 


4096 

i 

32769 

l 


c 

C40ort> 

occcoo 

8340000 

GOCOOO 

4395 

037401 

017,5 31 

143233 

1 7421 7 

4C9S 

0374 C 1 

0C7 177 

037401 

174217 

4096 

037400 

177532 

115160 

1 74 217 

3 192 

036020 

C 22 2 1 2 

1 72764 

171167 

•31 32 

036020 

00 OS 31 

G342 73 

1 71167 

31 92 

03601 7 

15 3207 

014117 

1 71 166 

1 2 2 fl 8 

C3352C 

C 34 C 56 

C31C73 

166777 

1 2288 

0 3 35 ? n 

CC 1 314 

111036 

166777 

12286 

033517 

157312 

16614 7 

166776 

153 8 9 

020372 

020225 

0 3 54 1 3 

1 65535 

3 5 38 9 

C 303 7 1 

154503 

175326 

165535 

IE 184 

030271 

1 27777 

CC4G3r 

165535 

2C480 

024531 

11 544? 

147551 

165261 

s G 9 a □ 

024531 

041470 

06342? 

165261 

2C«PC 

024531 

C 14 23 2 

07411? 

185261 

24576 

020311 

151353 

024561 

1 66003 

29576 

0203 1 1 

C 6621 4 

D74 166 

1660G3 

24576 

C2031 1 

042222 

16714? 

166003 

23572 

013661 

03 7430 

023745 

167502 

28672 

Cl 3660 

126277 

C 3 2 0 2 0 

i&rsn 2 

28572 

01 36S0 

105604 

163C73 

167502 

32765 

CC71<7 6 

G 35 51 7 

072200 

172111 

3?75 : 

9G7173 

1 4 C 3 2 6 

Oil 357 

172111 

3 2 7 5 : 

0071 7 9 

124551 

027731 

172111 

-4256 

oi3S5 7 

1C2727 

C50361 

1675C3 

-4Q9' 

01 3557 

173761 

1C40S7 

1 67 SO 3 

-4C5S 

013557 

17?76 I 

104057 

167EG3 

-9192 

C70 3 1 n 

047711 

151327 
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02031'' 

12336? 
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1 51 5 SI 
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024530 
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-16389 

030370 

174024 
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030371 

027547 
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030371 
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033517 

102546 
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033517 
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GO 7 55 5 

166776 

-24576 
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171166 

-24576 
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171166 
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-28672 
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157731 
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062092 
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01526.1 

C03 c 37 
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021650 
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023765 

CS&5GS 

021777 
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C5S677 
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CSS&61 
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0 to 32,768 pulses and return to 0 again. 
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1 3383 

033520 

056057 

06238- 

1 66777 

027423 

116668 
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1 56777 

027423 
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1 56755 

13288 

03 35 20 
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166776 

140163 
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114776 
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165535 
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121403 
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030372 

00 1 2 1 1 

0755 37 

1 6S535 

070112 
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067440 

1 57546 

16389 
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1 54503 

1 74 777 

165535 
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CE744C 
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024531 

142677 

143646 
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055013 

076473 
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0105=4 

30480 

074531 

066726 

002265 

1 6S261 

05501 3 

076473 
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042401 

3048D 

024 5 3 1 

041 470 

063031 

165261 

003530 

146055 

02576= 
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024567 
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117007 

3957b 
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1 12206 

147232 

166002 

172531 

024567 

031463 

033664 

0265 72 

24576 

030311 

066214 

073671 

1 66002 

1 31706 

□24144 

031«63 

033664 

02657" 

28672 

013661 

C 40222 

175246 

167502 

077256 

051726 

034414 

031327 
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33672 

01366-0 

147172 

062625 

1 675*02 

077256 

051726 

034414 

05 1 n 71 

1 71278 

28672 

01 3660 

126377 

031617 

167502 

G47514 

154136 

034414 

CP 1" 7 1 

171275 

327 63 

0G7174 

05 1271 

011163 

172111 

021154 

010755 

□36472 
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1 252P5 

32768 

0071 73 

154103 

165332 

17211 1 

021154 

010755 

03647? 

121110 

1S260? 

32768 

0071 73 

140326 

011272 

172111 
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074324 
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1526 0? 

-9036 

013557 

123521 

073700 

167503 
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170616 

034415 
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1 72 32 " 

-9C96 

01366C 

C 1 4 5 53 

050077 

167503 
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170616 

0344 14 

IS222C 

032033 

-9096 

013660 

014553 

050077 

167302 

111005 

- 064750 

034414 

131 42S 

1 45075 

-8192 

670310 

073703 

11277 1 

1660C3 

037402 

043077 

031464 

005314 

010441 

-8192 

070310 

157053 

□65000 

1 66003 

037402 

043077 

031463 

144472 

1 3546" 

-8192 

020310 

15 7353 

055000 

1 66002 

1 54232 

155704 

031463 

120900 

1 51 31 6 

-12386 

07453T 

064-304 

C56265 

165261 

062174 

115625 

02576= 

CC2=6C 

192767 

-13293 

1174530 

140260 

012345 

1 6S261 

0621 74 

11S62S 

0257-9 

131?77 

1 1 6203 

-12288 

0745 30 

1 40260 

012345 

165261 

006221 

060264 

025765 

1 1 4 " 4 2 

0=7662" 

-16339 

030379 

02 05 32 

003627 

165 S3 5 

035406 

142665 

021651 

102=22 

0=1 1 == 

-16334 

030371 

06 4254 

130311 

1 65 53 S 

□35406 

1 4266 S 

021651 

021=37 

144137 

-16364 

030371 

064254 

130311 

165534 

171664 

123674 

02165" 
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-20432 

033517 
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0144 SO 

16S776 

1 36227 

146406 
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-2P48C 
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124550 

023216 

16617 76 

136227 
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1 0346S 
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015261* 
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-29576 

03601 7 
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135232 

171166 

073215 

073755 

010571 
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15235= 

-29576 

0 "601 7 

141624 

0021 35 

1 71166 

□73215 

0737SS 
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064322 

07027" 

*29576 

036017 

141624 

0021 3S 

171166 

051632 

140535 

010571 

046=77 

93534 " 

-28672 

037400 

156733 

034610 

174217 

000624 

013211 

004 IS" 

017223 

0 "60 P" 

-28672 

037400 

167272 

174254 

1 74217 

□00624 

013211 

004157 

120221 

=157555 

-28672 

0374 CO 

167272 

1 74259 

17=7216 

170264 

16211 1 

DC4157 

110657 

05066= 

-32768 

037777 

177777 

037523 

177777 
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037147 

000009 

1 20"0 1 

0 3 7 65" 

-32768 

037777 

1 77777 

1 37525 

177.777 
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037147 
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02765" 

-32768 

037777 

1 77777 

1 37525 

17 7*777 
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037145 

DOOOOO 
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03765" 

-32769 

040000 

ooaooo 
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000000 

DOOOOO 

037144 

00000=7 

000000 

=13 7 653 

Three 

axis Az, 

i Ax, Ay 

incrementally - 0 

to 32, 

768 pulses 

and return to 0 

again. 
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Single Axis - Runs 1-5 
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Oscillatory Runs-Ten Cycles 0 to 1024 
Pulses and Return to 0 Again 
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Single Axis - Runs 6-10 
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C;7 7 70 

C*. 02 4 7 

163131 

COCQOC 

-512 

CZ 777C 

CGCC02 

163404 

OOCOCC 

-1G24 

owoior 

COCGGP 

C3520C 

000000 

-612 

ruccrr 

ccccor 

C3C 200 

COCOCO 

-1324 

"43Grr 

C37CG1 

036200 

000330 


L 

r4pnrr 

rcr:c° 

C 7 6 20 0 

DOODDC 

512 

r u o c 30 

oc'ioc 


COOOOO 

1324 

On rr~ 

cocoa' 

CJ52C0 

u C 3 30 0 

512 

C 7 7 7 76 

CCTC c° 

163273 

000000 

102 4 

O’ 777 1 ' 

OCDC57 

16300 3 

COCCGC 

512 

C 3 7 766 

0CC15S 

117272 

177400 

1 


CC~525 

1 G'>4 □ 5 

1 77030 

12 

C 7 77f*r 

CO" 15 = 

10 7 23f 

17 7400 

-1Z7** 

■•>37770 

CGC - 5"* 

1 5^5 S' 

GG03C0 

-■ 12 

C 7 777c 

LCO'O 7 • 

162 1-4 

CCC3GC 

- 1 “24 


GO 233C 

G357 S4 

003DCC 

- : ! I 

r l*i"~ ~r 

: ccro r 

C35754 

COCQOC 

- n?-4 


w Z T"" 

036754 

0 oooon 


z 

Z 4 


C05754 

g ccoon 

riz 

njjrv"" 


7 5*4 

G DCOOC 

1 "4 

*i ~ “ r *" 

3’*' 

. 3 = 7 54 

cooair 

512 

3 3 7 7 76 


1 o3Gc7 

aooaoo 

10 2 4 

p 3777f 

GCC o S 2 

1 6 7 52 5 

aococo 

1 2 

f 3 7 7 4 r 

t C C 1 5 6 

1C7CG7 

1774 jr 

l'?4 

7 7 7 76** 

00C526 

1 C 7 1 75 

1 77GOG 

-- 13 

f ' 77T C 

'•»' 1 C 

«. - J 

U^77 c 

17 7400 

-1C24 

*“ 3 77 7“ 

.3"_j 7 

1 .74 16 

-uOOIG 

-512 

T 7 7 7 7f 


I S ? 7 u S 

CLCGCC 

-1C24 

(3«i nr or 

33 3030 

03 5 5 4 r 

ooooor 

-5 12 

rm-rcr 

t or:cr 

C 35 54 C 

cccoor 

- 1 : zn 

c ^ g r “t r 

‘03.3- 

4- Z c E 4? 

3GG1C0 


coococ 

C34526 

DGOCOC 

ooooco 

C466rn 

ccococ 

l.34626 

ocoocr 

GOG°DG 

C466=T 

GDDGCC 

034=26 

00000" 

000 n 0G 

04660' 

coococ 

C34S26 

G00377 

176=25 

1 225 3 7 

oooeco 

034S26 

000777 

165 7 S3 

C 36 52 1 

021251 

1 606 20 

000777 

165’53 

036521 

G52521 

04=1 73 

OGO 7 77 

165753 

C36S21 

G 2 1 2 5 1 

160540 

000777 

1 65 7 5 3 

C 36 571 

COCOCG 

C33C S3 

000777 

165 7 53 

C36521 

aoocoG 

J336S3 

□ GO 3 7 7 

176=25 

1 23 41 7 

GOOOOG 

333653 

oaaooo 

OOO n GO 

P 50 33 3 

CPOCCC 

C33 6 53 

ccocor 

000PC1 

C5C3 3 7 

OOOCQO 

033653 

00003" 

OOOnOQ 

■ISC 33 3 


□03C0C 

O 7 3 6 5 3 

00000' 

000°03 

"50 33 7 

-• <\r' rs f' 

Li U U< _■ k# 

0 3 7 6 5 3 

OCCCO' 

ooo°co 

" 50 33 7 

LQCCGC 

033653 

ooccor 

GQC°C C 

rFC3i~ 

OOOCQO 

G33653 

000377 

176=25 

1 24277 

CGCCCO 

633653 

CC0777 

165753 

C4025 1 

02X251 

1=7753 

COO 7 7 7 

1 65'53 

040251 

CS2521 

G44326 

DC0777 

165753 

=00251 

62 125 1 

1=7673 

CCC77 7 

1 6 5*>5 3 

=4C251 

GOGCOO 

033002 

000777 

IS5 7 53 

=402=1 

CCCCQC 

C33CC2 

DCC377 

176=25 

12S16 7 

OOQCOG 

033002 

QCCOC' 

OOG'GO 

"5237' 

CDCCOC 

03 = 00,? 

oCCOOO 

ODC r 06 

C52G7= 

GCG1CC 

33 3C0'2 

GC 001' 

OOO'OO 

05207- 


cccrcc 

0I3C02 

OODDOP 

ODOnCO 

C=2G7" 

COOC0G 

0 3 7 3C 2 

ocooor 

OOO'OQ 

-5207" 

QOCOCG 

03=002 

00301° 

003°0C 

C 5237' 

oooocg 

333002 

000377 

176=25 

1 2602= 

CODGro 

C33002 

CC0777 

165753 

0420-14 

D212S1 

157103 

000777 

165753 

0420°4 

05ZS23 

043461 

00C777 

1 65 7 5 3 

r 42004 

02 1 26 1 

1 = 7C 2 3 

CC0777 

165753 

' 4 2C °4 

oaocoo 

0321 3J 

COO 7 7 7 

165753 

" 4 20°4 

CCGCCu 

C 7 21 3. 

0CC277 

176=25 

1 7 67 C 7 

oooccc 

3371 3. 

OCOOOO 

oconao 

P 53 575 

COOCCL 

3321 36 

GCOCOn 

dcg'go 

05362= 

COOnCi. 

J 3 7 1 30 

OQOG3° 

ooanoo 

r 63 62 = 


aooooc 

3371 33 

OGCOir 

GOO'CG 

C 5367 ; 

CCO'C. 

33713. 

030C1" 

ooo'og 

“ 53 I 7 ^ 

6 330'C 

037130 

OC GCC° 

O°0'OC 

0 =3 52 = 

ooococ 

C 3 2 1 30 

OGO 3 7 7 

1 7 5=2 5 

1 2755 7 

CGCCOC 

6321 3G 

COG 7 77 

165753 

04353' 

C 2 1 2 5 1 

is r : 7 3 

GCC77 7 

155753 

= 435 37 

G 525 ? 1 

54751o 

DC 0 7 7 7 

165753 

043537 

C?l?= 1 

1==1 S3 

LCP777 

1=5-53 

043537 

occ'r. 

j31 25o 

0C07 77 

155753 

'43537 

COCCGC 

63125b 

GCC277 

176=25 

1 7 04 3 7 

OOCOCC 

0 7 1 2 55 

00000° 

Q3G1CG 

,155355 

LGC1CC 

C 3 1 256 

c c n c nn 

crG°cc 

C=5356 

oco-oo 

031 256 

c:gio° 

OnOiCI 

' 5=35= 


Three Axis - Runs 5-8 
TABLE 5.3 (cont.) 


120 



0 

naoonn 

OOOOOP 

035548 - 

&03000 

000000 

031 256 

OOOCOO 

ooonoo 

055356 

512 

04ccrr 

ccccgo 

03 4 54 0 

000000 

000000 

t31 256 

ococcr 

occ'cc 

r'sssp 

10211 

040000 

cocoa" 

035540 

000330 

000003 

3312S6 

OCCOG n 

aoa r ’ao 

<'45355 

512 

037776 

3CC332 

152572 

Qoaaoa 

QOOCGC 

031 256 

0CC 3 77 

175425 

1 31 306 

10211 

037770 

000052 

IS 22 56 

ooaaoo 

oooooc 

□31256 

000777 

165753 

C 45256 

S12 

R3776P 

OCCIBU 

106522 

1T750C 

021251 

155364 

DC0777 

1657S3 

C4S256 

1024 

□3774" 

COO 5 25 

10 1 5 77 

1 77000 

052521 

041750 

000777 

165753 

045255 

-512 

0 7 7 7 6 P 

CG315 C 

1 C 6 50 7 

1774CC 

021251 

156304 

CC0777 

15S7S3 

C4EZ56 

-1024 

r 7 7 7 7C 

C3CC5 " 

1521 36 

OCCOOO 

OOOOOC 

03"403 

000777 

155753 

T 4 5255 

-512 

C? 7776 

CCC23 ? 

1C 2 44 6 

i* C 3 C C C 

CCCCCC 

031-403 

CC0377 

I7£*25 

132167 

-1024 

ctoccr 

coccr”: 

Ci ir 314 

CODOOO 

COCCOO 

030403 

GOCOCC 

goc'gc 

r 6 7 1 06 

-512 

o 4 o o o c 

oocoo n 

C35314 

000000 

OC0"0O 

030403 

OOOCOO 

OOD ri CO 

0S7106 

-1024 

04pppp 

ccr:c r 

C35314 

oocaoc 

CODCOC 

03T403 

ccocc r 

OCC'OC 

05 71 "6 


i_ zc r zz~ 035314 aoooco codopd C 32433 aoocon a^o^QC rsii"'- 


412 

04007" 

LC"CC" 

035314 

OOCOOP 

oojrcc 

G3"433 

□0000" 

□"□"23 

rs?i"4 

102-* 

cecrcr 

icc :cr 

L:=31 4 

cococr 

cccrrc 

C3 C4 0 3 

cccccc 

CCC"CC 

r 5 71 "6 

512 

C2 7770 

CCCCC7 

162336 

OOCODO 

cccccc 

G3"4D3 

GCQ37 7 

n&'as 

1 33C 3 C 

1024 

03 777T 

OC0O57 

161777 

000000 

DOQOCG 

377422 

003777 

155753 

C473C1 

512 

037766 

crci55 

1 G 5 23 5 

1774CC 

□21251 

154512 

0CQ777 

165753 

C 4 70 0" 

1024 

037760 

CC C 5 2 C 

121406 

177000 

052521 

041 101 

QC0777 

165753 

3 4 7"<” 

- 51 2 

r 3 7 7 6 b 

CCC15= 

106222 

177400 

C2 1 ? c 1 

154432 

CCC77 7 

165752 

C 4 73-C " 

- 1 C24 

037770 

CC0C52 

161654 

ooaoao 

0 030 CC 

3275 33 

000777 

165753 

C470O" 

-512 

037776 

CC0CC2 

1S22C3- - 

CODOOO 

CCCCCCi 

02 7530 

DCG377 

17&<= 25 

1 7 37 3 7 

-1024 

0400110 

000300 

0350 61 

000000 

OOOCOO 

□ 27 5 30 

OOOCOO 

QP073D 

0 SO 64 ' 

-512 

C40C0" 

CCCC30 

C35C61 

oooooc 

COCCOO 

C 2 7 5 3 0 

CCODOO 

CGGPGO 

C&06 4 - 

-1C24 

04 OCT" 

0CC3C" 

C3 C 361 

ooooon 

OOOOOC 

027533 

000000 

330 n 0a 

7 5C 64" 


Three Axis - Runs 9 and 10 
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3F0=»* I S JUMTR 

CYC LF 0 , ' n CCmtEO BY 1201 00B7D "N " > 7 FE8 70 AT 10:*H7I52- 


M» is ^ rV OB *? A M 


BTf- ar.L 

ussc : 

coon l > 

ocrrm 5 i 

0 A I A ( 0 1 

OQOG42S 

Blank con«Q»<2l 

OOOOO 3 









m l 


r^’75 

irn 

cm 


crrcrF 

1 ZL 

crci 


coon c 

1S0C 

COO! 


COCl’2 

1 5 e c 

com 


nor i3» 

1F4C 

rri 



1 TIC 

cnn 


COC „ T 7 

1 ELL 

cm 


000777 

1 a 9L 

0CC1 


OOCIM 

2CCC 

0C01 


000*73 

?rci 

rcc i 


racers 

7 90L 

□ GDI 


Lr0075 

4 39 L 

C0C1 


COG 4 71 

300L 

0001 


OQC 3 C 5 

3025 

com 


0004C3 

3 1 3C 

"c: i 


pr?t*o s 

376C 

ocn 


CCC4 54 

7 3ZG 

ecu 


000463 

3 36 C 

C001 


noCS02 

35CE 

oom 


000473 

3* 0L 

non 


rr p r 74 

4356 

cm 


CCC7 t K 

446C 

0001 


000756 

4 6 1 C 

0001 


nor7R3 

46 5 5 

com 


COOT 7 c 

9 71G 

rrr i 


rr^rMe 

4°Ot 

nnn] 


C0C7 nR 

499L 

cm 


000702 

b ro l 

cool 


000746 

59TL 

onoi 


t)0nT74 

S0GL 

rr' l 


m fi j ror 

P PCL 

core 


CO.OC2S 

4 1 CO F 

coco 


ccrc?6 

AJOCC 

0001 


PG0241 

9GL 

anno' 


000027 

° reef 

r * 


tt:’ 

eiriF 

cror 


rrcc ’2 

11 2 cor 

crco 

I 

000004 

oil 

DOCC 

I 

000007 

C 1 2 

uono 

i 

onori? 

C 13 

0000 

I 

cporni 

0* 

u C C o 

1 

000002 

:i 

C0U0 

I 

Q00C03 

n: 

aaoo 

I 

000015 

IN 

oanj 

i 

annm 

It 

0002 

I 

rc n ?r ^ 

J 

aoac 

I 

000024 

K 

GOOD 

I 

000023 

L 8 

0000 

I 

0000 ’2 

LIM 

0000 

i 

ononis 

LIM1 

mor 

i 

# *' i 2'M 7 

Li v i 

oaoo 

I 

303 u2C 

tIM J 

0000 

I 

0 00000 

OUT 

0002 

I 

nooa n a 

03 

0002 

i 

mnm 

7 


Main Program 


TABLE 5/4 

Program Listing 



Coll ’ 

3 * 

0 - 

-C5“i T( 1P1 } 


rc 11 7 

10 - 

* i" i 

-ORf* T (/) 


<-c 1 1 1 

1 1 - 


Cl T 1 ,. 


"Cll' 

l 7 * 


:\ :1 


Ci! 1 1 * 

10 * 

i’ 

Of 4 0 t T i< • s 2 C C I CX*C.YiOI 


00 1 7 3 

I 4 * 


IF I CX.CS.-l ) 5TC ' 7 


rci- c 

l c * 


»T I T E ( C,J T 1 r P r C ) Z X . TY ... ' 


n; i t 7 

13* 


7 £3^ ( 1 'J « a 2 C 0 l P-ifLO' li O *’ -7 


''ciur 

17* 


w rtt e t gut , ? ?rr. ) FR»Lii.i*Li**-»L:r’ 


f’Cl 43 

1 a * 


i- ao t in, ~ 1 no 1 '"ii 


r D i r u 

1 c * 


?E4 cc :> • e 1 cc i c 


0013? 

20 « 


31 401 IN. 5100 c 1 3 


02 1 7 F 

21 * 


DC 7CC LT4-1>I.I1’3 


HQ17 7 

2 7 * 


L3-C 


"Cl 7 <* 



2-1 


f 'Cl 7 - 

_ -t * 


VO ITE ( CU' . 7 OOC 1 l2»:11>C1?.013 


93 21 « 

2 "* 


CO 3E0 K - 1 « U IM 2 


rm 7 

7- * 


i-:?4 *? 


OOZ?" 

27* 


IFlCX.E3.01 GOTO 13= 

„ 

CO 22? 

2 « * 


0 0 ICC Jzl.LTM 


9C2 7 5 

2 7 * 


C 4 LL IT£?{C12iC13) 


n 0 >7 r 

" * 


I c l (LJ V 1 .t’~. ! ) .4 NO. 1 1 CO I* ■*f).F u .E)l 

0 0 T c 0 " 

I' ^ -I ■» " 

01 * 

■ 

1 'C'CCtj* 331,4.5.0) r C '0 13. 


"C 2 -- 

22 * 


L 2 : J * 0 


rrjT* 

32 * 

°c 

*Rir EicuT .ircco Lu»ciii vio.cio 


00232 

34* 

1 51 

CONTIS J£ 


n C 2 = 4 

35 * 

15 c 

I F (CY.EG.n CCTC 2°S 


■'0C3P 

Jk * 


:c oc? j=i*Lim 



17 * 


C 4 Ll I XE= l C 1 7 iC 1 1 1 


*023? 

Z** 


!r((Li f 11. r “..l) . -”C. I '0 0 ( <• 7 • 1 . .0) ) 

JCTJ J?' 

r C 2 F« 

z\* 


IF (. V CDU»F~ ) .Nt .0 :c~: oco 


00233 

U" * 


L2 -J* 7 


9C2S7 

C! -r 

nr 

UR I T E ( CUT . a CCO 1 Lc • C 1 1 . ClO . 2 1 1 


i** ’n' 

«■>* 


C C NT I SCIf 


* 0:1 * 

43 * 

7 - ' 

iMor.EC.o 3cic 3?r 


* 0:1 ? 

Zf £l * 


jo orr ^=i » l i v i 


"C 31 3 

45* 


C4LL IT£?tCil.Cl t ) 


"C Z 1 5 

45* 


IFdLm.cO.l ) • 4 ,N C . I'-O 0 1 < 1 7 5 1 . "j .on 

30 T 0 ? 7r i 

*0 i “* 

47* 


7 F (reel J tFR > .NE 1 f v tc i-o 


*" 1 7 7 

4 3 * 


L a :j*7 


f*^'77 

4 

"* j “ 

-7 it e ( cut ( ° r cc 1 Lc « r 1 : • cir . ;i : 


"Ola? 

j" * 

35"' 

: omi,' o: 


3 L u 

51 » 

5 3 5 

Of Til L£ 


•'Z !« = 

= 


» “ _ 1 


''033 7 

£ 3 * 


re eEC Kr 1 » L If 2 


3 c 7 

5 4 * 


L n -5 * 7 


'■"in 

5 5 * 


IP (102.-3.0 :ci r O' - ' 


^ ^ 3 ^ ~ 

z £ * 


03 !i:T J:1 iLI' , 1 


*c:' n 

5? « 


0 4LL IT£“ IC11 iCl tl 


n C j31 

3 4* 


i-<(limi.fc.i) . 4 so. i-,:;ur 5 'i. c ;.;n 

GOTO 77 0 

r G3 :7 

3 c « 


IF tf CC ( J (Pi 1 .NE . 0 ) GOTO *.cc 


7 033 c 

t>" « 


L a =J * 7 


r;»'r 

C I* 

7 — «*■ 

7>9iT*(rj T i 3 rcci 1 


'’04* = 


45" 

Z C v T I *J J z 





Main Program (cont.) 





TABLE 5.4 (cont.) 
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rc-.~7 

'2 ♦ 

: 3 " : P ( : y . i ^ r > : c r c 4 ss 

'’C.i 1 

34 « 

Cf j:l .L Il'l 

~C 41 4 

6E * 

CALL I r £P (Cl 2iC 1 1 ) 

PC4 I S 

pfi * 

!f(ai' , i.M.n .a no. woji^.o 

r £ '!.i 7 

67- 

; f i.'::ij(F"i.fE.oi rou 

"C-t ■>! 

3»* 

L3-J-7 

vz^->? 

6 3* 

“SC h £ ?IT3[CU T t'CCC) LctCUrC12.Ci; 

r^“4 1 

7C* 

SC E CONTINUE 

cc-* 1 *-? 

71 - 

439 IrtCX.C3.Cl 5 0 TO 65C 

Pf 44 e 

7?* 

cc t cr j ; 1 • l ip 1 

r ':‘tr 

73* 

CALL ITC71 Cl 2.C13 > 

rC4 « i 

74« 

:=■ 1 1 l i;- 1 . 1 j • a nc. (^su^ri.Fj.m 

p "> 4 = 7 

7^* 

: r t u cct j. ? 3i .ki.z ) oo. 

r''ij:c 

7 r * 

t -:jr 

""4 

77. 

■ :t~ t :-j* . =icOu i l 3.c 1 i.c 1 :» ~ 13 

rr4-> = 

7 ' * 

-r- CONTINUE 

n C4’7 

7 3« 

3“ CONTINUE 

rc;°i 

et* * 

* c I TE 1 CU T . ^206 1 

p ~c - 3 

'.1 * 


"C c "5 

4 r . 

in. ; ; 

n -j ; n ; 

4 3 * 

: ^ 

Z - 

c " u -il a ti o : f.o oi&G^o ^Ti:'' • 


Main Program (cont.) 
TABLE 5.4 (cont.) 
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3F03.I- Mfl T 

'■yet- -: r C^^U'" *V 1? 1 CC<70 PN ■>7 FEd 7C St 1i_:4P:£f. 


•t ■ > ' C UT ! K F I T E? 


y =rm cco3i2 


'tTfutrE i i cfo: cf' c m rnL3^:. o a x a t o » ccocti: hssk cofmcmd gcgoci 


FyT r 3NSL •’ r Fr ! »:Nr‘' S (ELOCrti M“£) 

OOP: 'CALF 

0004 NEREIS 


t STORAGE A5SIGNHENT IdLOCKt T7PE > RELATIVE LOCATKN, NAME! 


ro 

0000 


Q0CC35 

9000F 

0000 

r 

D 00010 

IEC 

COCO 

I 

GD0CC5 

IES 

OOOO 

r 

DC r C 0 6 

IFC 

OOOO 

I 

000023 

IF 5 

cn 

rcoo 

i 

000075 

TfC 

DPCO 

i 

00007? 

IGS 

crco 

I 

C0PC24 

IRC 

0000 

I 

rorc’i 

IRS 

anno 

t 

nnnoi? 

I JC 

i 

oooo 

i 

000077 

I JS 

0000 

r 

000033 

IK C 

0 00 0 

I 

000030 

IKS 

oooo 

1 

OOOO 34 

ILC 

oooo 

I 

OODO 31 

ics 


0000 


000040 

I NJPS 

0000 

i 

□00011 

IRC 

0 00 0 

I 

000006 

IRS 

oooo 

I 

000017 

iyc 

oooo 

I 

000007 

175 


OQ 04 

i 

00 QQQ 2 

J 

0000 

i 

DOOG 1 & 

JEC 

LOCO 

I 

QD0C13 

JES 

OOQJO 

1 

000017 

JflC 

oooo 

I 

000014 

JF5 


fioao 

i 

000070 

JtC 

GCDfi 

l 

DODO 15 

JVS 

DO 00 

1 

DD 0002 

OR 

oooo 

I 

oono n 3 

HN 

□ood 

i 

oocca i 

r»t 


0000 

i 

000004 

Bfi 

4 030$ 

t 

000000 

OUT 

0007 

I 

000000 

PR 

0002 

I 

000001 

7 

• 





Subroutine ITER (performs one iteration on two direction cosines) 


TABLE 5.4 (cont.) 



-omi 

! « 

5UB' J CUTIK.C ITESIC.C'r) 

031G3 

1 « 

PjTE 311 PlfZfJiCtCYiG 11 

'-cunti 

^ • 

2 IHf A ' row C ( 31 .CYt 3 1 

001 n S 

* * 

CC^'O*- ?7,7,j 

"o i r-f 

c * 

OiTt i/.'K/f's, />• 

00113 

J • 

3AT« CJT/5/ 

rci l c 

7* 

-rcr 1 cof « 7{ r ip, & (><x , os 1 1 

"1U C 

3* 

i-:o:ct 3 1 

00117 

!♦ 

T4S-C( 2 1 

cci?r 


tYS-CI 1 1 

001?1 

1 1 * 

:ec;cy ( 31 

CO 122 

12- 

1KC;CY(21 

00123 

13. 

IYCrCV( 1 ) 

001 7 U 


JES: [CS/rC 

00125 

Tt 

JSS=I 7S/V0 

0C1 ?F 

1 - * 

JYE: I Y s / ^ r 

OQl’7 

'l 7. 

jec-i':/ u ; 

C013C 

\ fc * 

J6c~i c c/ M r 

00131 

1 ^ 

JYC:I v C/“0 

C01 72 

2 t * 
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Subroutine SCALE 

(scales one direction cosine by operating 
on each of the three integer parts) 
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APPENDIX A 

SMOOTH PULSE SEQUENCES 
( SRRC-RR-68-48 ) 



SMOOTH PULSE SEQUENCES 


A. J. Lincoln, M. Cohn and S. Even 
Sperry Rand Research Center, Sudbury, Massachusetts 

INTRODUCTION 

The necessity sometimes arises in digital systems to generate ex- 
ternal pulse rates whose- frequency is related to the internal clock frequency 

by some proper rational fraction while still maintaining synchronism with the 

* 

internal clock source. In many such cases, it is a requirement that pulses 
of the generated frequency be as uniformly spaced as possible within the con- 
straint of synchronism. We define such pulse sequences as smooth pulse se- 
quences in contrast with uniform sequences in which each pulse is separated 
from its predecessor by a fixed interval. . The smooth sequences possess a 
number of interesting properties, both in their structure and in the ways they 
can be generated. Some of these properties will be treated in the following 
discussion. 

Although it has not, to the authors' knowledge, been previously 
identified as such, a method has been in use for some time for generating 
smooth sequences. This is the constant multiplier used in digital differen- 
tial analyzers CDDA's). Another technique often used for forming synchronous 
frequency multiples uses the "rate multiplier," and involves detecting the 
one-to-zero transitions of the stages of a binary counter chain. The rate 
multiplier technique, although capable of forming relatively smooth sequences, 
does not form a true smooth sequence, as can be easily demonstrated. 
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Other methods for generating smooth sequences can be described. 

Our interest here, however, is to provide a formal foundation for the subjec 
and to prove several important properties of smooth pulse sequences. 

PROPERTIES OF SMOOTH SEQUENCES 

Choosing an arbitrary reference point in time, and assuming, for 

convenience, a unit-period clock, we define two types of pulse sequences: 

Uniform Sequences and Smooth Sequences. Each can be described by a doubly- 

t h 

infinite sequence of real numbers, the i L number giving the arrival time, 

"L h 

as counted by the clock, of the 1 pulse. Since our reference-time zero 
and the position of the index zero can be chosen arbitrarily and independent 
iy. we wish to consider equivalent those sequences which differ only by a 
shift in time or by reindexing. In other words, only the "pattern" of a 
sequence is important. Accordingly, we first make the following definition: 

Definition 1 : 

Two arrival-time sequences {a^} and [b x } are equivalent 
if for some integers t and t , 

a. = b . , + 1 . 

1 1+T ' 

A uniform pulse sequence is an ideally smooth sequence with rate 
p/q times that of the clock. It consists of p uniformly spaced pulses 
during every q clock periods. 

Definition 2 : 

A Uniform Pulse Sequence of rate p/q has pulses at 

times {u } , where u. - i i . 
l i p 



Unless q/p is an integer, the pulses of the uniform sequence do 
not all coincide in time tv 1 1 It clock pulses, so the problem of synchronous 
pulse rate generation cannot generally be solved by <i uniform sequence. The 
concept, though, is crucial in defining and character! zj ng smooth sequences, 
which, heunstically, are the best synchronous approxinat ion to uniform 
sequences. 

Definition 3 : 

A Smooth Pulse Sequence of rate p/q has pulsus rmJy at 

integral clock times {s^} » - ^ i - ; moieover, these 

pulses can be put into one-to-one correspondence witn the 

pulses of the uniform sequence in such a way that 

raax(s.-u^) - mints. -u ) < 1. 
i i 

To construct a concrete example of a smooth pulse sequence for tne 
fraction 5/8, we place pulses at those clock times which coincide with or 
immediately follow pulses of the uniform sequence. 

Example 1 


Time 

r— l 

CU 

1 

0 

1 

2 3a 

5 

6 7 

8 

Clock 

^ 

_A_ 

__/V 

AAA 

_A 

A A 

_A_ 

Smooth 

A 



A A 


A 


Uniform 

/ 

/ 

A.. 

t 


/ ^ 

A A 


/ 

A 

1 


0 /»■) 

+ 

8/5 

-®|*- 8/5 8/5 

-X. 

i 

8/5 0/5 

i 


Here the arrival times of the smooth sequence are >..., -1, 0, 2 , 4, 5, 7, 
8, ...} . If the pulse pattern between, say, time 0 through time 7 is 



repeated indefinitely, an infinite, smooth sequence of rate 5/8 results. 
Since this sequence was constructed to be periodic, it suffices to observe 
the smoothness of a single period. Under the natural correspondence (shown 
by arrows), the deviations between the smooth and the uniform arrival times 
range between 0 and 4/5 , satisfying Definition 3. 

The last condition in Definition 3 bounds the range of time devia- 
tions between pulses of a smooth sequence and their correspondents in a 
uniform sequence. We next observe that the bound can be sharpened and that 
within the equivalence relation the deviations themselves can be bounded. 

Since the Sj axe integers, and the u i are integral multiples of 
q/p , the differences (s^-u^ between corresponding pulse times must be 
integral multiples of 1/p , hence 

Lemma 1 : A smooth sequence satisfies the inequality 

max(s.-u.) - mints. -u ) ^ . 

,ii i i i p 

Now suppose that i is an index where the difference (s^u.) 
achieves its minimum;^ that is, 


mints -u.) = (s. -u. ) 

^11 1 0 1q 

By an integral shift in time and a reindexing, we form the new sequence {s*} , 
where 



To be rigorous here we should note that Lemma 1 and its preceding discussion 
imply that the differences (s -u.) take on at most p distinct values, and 
that these values are bounded iielow by Sj-uj -(p-l)/p, where j is an index 
chosen arbitrarily. Therefore, the sequence J has a greatest lower bound, which 
is s. -u. for some i. 


s. -u. 
*0 io 


0 




The minimum deviation of [s.'} from the uniform sequence is given by 


minCs^-u.) - minis.. . -s. -u.) 


i 1 


i + in io i 


minis. ,. -u -u. -s. +u ) 
. i+l 0 l io io io 


minis... -u. , ) - is. -u ) 

. i+i 0 i+i 0 io io 


But by definition of ig this quantity vanishes, so that all deviations of 
{s^} from the uniform must be nonnegative. Invoking Lemma 1, we can now 
state : 


Lemma 2 : 


Within equivalence, every smooth sequence satisfies 


0 ^ is -u . ) ^ 

i l p 


. eo ^ j[ 5 o 


The sequence in Example 1 was designed to satisfy this bound without an 
equivalence transformation. 

We are now able to prove the fundamental theorem of smooth sequences. 


Theorem 1 : A smooth sequence of rate p/q is unique to within 

equivalence. 

Proof : Given any two smooth sequences of rate p/q, consider their equivalent 

X L 

forms complying with Lemma 2. The i pulses of both smooth sequences must coincide 
with clock pulses, and must coincide with, or lag by less than a clock period, 
the same correspondent, u , in the uniform sequence of rate p/q , There- 
fore, the pulses in the two smooth sequences must coincide with the same 

clock pulse, hence with each other. QED 

This is an extremely strong result, following, as it does, only from 
the definitions ol uniform and smooth sequences and the notion of equivalence. 



So far we have used arrival times to describe pulse sequences. Tt 
is convenient now to introduce a dual description, the cumulative count. If 
A(t) is the number of pulses in the sequence [a^l which have arrived by 
time t , we define the "count sequence" (aKt)} , where 

aCt) = A(t) - A(a ) 
o 

Thus aCt) is an integer-valued function of time, satisfying 

Lemma 3 : If {a.} and [aCj)} are, respectively, the arrival 

time sequence and the count sequence for a synchronous 
pulse sequence, 

«tJ> = i for * J < a i+1 
a or( j ) ^ ^ < ^(j Hi 

As done in Lemma 3, we will always use the lowercase Greek equivalent 
of the arrival-time character for the corresponding count sequence. Also, 
by analogy with the notation {u } for the arrival times of pulses in the 
uniform sequence, we make the following definition. 

Definition 4 . A uniform sequence of rate p/q has a (rational) 
count sequence {u(j)} , where 

u(j) = j ^ , - «> < j £ » . 

We are now ready to prove the important relationship between the 
arrival time sequence and the count sequence. 
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Theorem 2 ( Duality Theorem ) : 


If ( s ^} anc * j )} describe a synchronous pulse 


sequence of rate £ , 

q 


max(a.-u.) - min(a.-u ) £ 
i 1 1 i 1 1 P 

if and only if 

max(u(j) - o?( j ) ) - min(u(j) - cKj)) £ 

J j q 

Proof: To prove the "if" part, observe that for any a. 

p m * nCu( j ) = p min (jj j - or(j)) £ -jj a ± - a(a.)) = ( a. - 

J j 


Also, since 

max(i 
P 


iax(u(j) -a(j) ) (a -13 -a(a.-l)) = (a. -l- q (i-1) ) , 

j p v qi i/ip / ‘ 

~p)' r f max(u( j ) ^ (a i“ u i ) • 

From these two results, 

V “n ) + « max(u(j) -«(j)) &max(a.-u.) , 


min(u(j) — o?( j ) ) < min(a.-u ) 


l l 


By subtraction and hypothesis, 

(- 1 ~ p ) + p ^ S maxCa.-u.) - min(a.-u.) 
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To prove the "only if" j|art, observe that for any clock time j , there is 
a pair of successive arrival times such that a £ 1 *' a. ,,-l , so that 

t 1 J i+l 

rv(j) - i . Therefore, 


q a i " 1 " q J " aCj) £ 
q ^ a i~ U i^ “ u (j) - 


£ (a. ,,-l) - i 
q i+l 


(a. ,, -u j.,) --^ + 1 
q i+l i+l q 


But then 


(l £ max( 

v q ' q , 


£ 

q 


a ,, -u ) 
l+l i+l 


mmCa . -u. ) 
i i 
i 


^ max(u(j) - cvC j ) ) 
J 

£ minCu(j) -cKj)) 

j 


t 


t 


again by subtraction and hypothesis, 


(l £ (j 2 — ^ ) = s max(u(j)-of(j)) - mm(u( j)-a'(j) 1 . 

q q p q j j 


QED 


An obvious consequence of Theorem 2 and Definition 3 is 

Corollary 2.1 : A synchronous pulse sequence described by the count 

sequence {<r(j)} is smooth if and onl\ if 

max(u(j) -g-(j))- minCu(j) -o'ij)) ^ 

j j q 

In the discussion which follows, several of the results have 
analogous statements and annloaous proofs in terms of arrival times as well 
as count sequences. In such cases, both statements, but only one proof, will 
be given. 


fir - 8 



Lemma 3 : 


-Li', for all i , a^ +c = aj +d , ihen for any integer k , 


Vkc = a i + k,J ■ 

L e [ri m a_3_ : If, for all j , crCj+d) = oKj ) + c , then for any integer k 

cvC j + kd ) = <v(j ' e ko. 


Proof s : 

Both proofs are simple 

inductions on k . 

Lemma -It 

If, for all i 

■ a i+cj = a i + d l 


and 

a i+c 2 = ' a i +d f! 


then 

c l^ c 2 d^ydrj 

Lemma 4 '. 

, If, for all j 

, cr(j +dj ) =ty(j ) + Cj 


and 

Q-(j +d i;) ) -Q’(j) tc 0 _ 


then 

C l /C 2 = ( V d ‘> * 

Proofs : 

Consider crCj +d^d 0 J ; 

applying Lemma 3' twice 


Oft jj + (3^ = Gf(j) + dg-rc, 

for all j , so that 

c i/ c 2 ~ d i/ d 2 * 

L emma 5 • The arrival time sequence fs ] for a smooth sequence of 

* 

rate ^7=^ , P and q relatively prime, satisfies 
s , = s +q for all i . 

i+p i n 

Lemma 5 ' : The count sequence foTj)} for a smooth sequenre of rate 

*77 = 7§ . P and q relatively prime, satisfies 
r(j-tq) = o(j ) + p for all j . 



Proofs : 


a(j+q) - r(j) = e(j+q) - u(j+q) - o( j ) + u(j ) + q ^7 

= (rCj+q) - u(j+q) - rrC j ) + u( j ) } + P ■ 

The difference on the left-hand side of this equation is clearly an integer, 
while the bracketed quantity on the right is less than unity (by Cor. 2.1). 
Therefore, the bracketed quantity must be zero, proving the lemma. 

We now define two new sequences, the first differences of the arrival 
times and the count sequence. The first, called the "gap" sequence, consists 
of the clock- time intervals between pulses. The second, called the "binary" 
sequence, is unity at those clock times when pulses are present, and is zero 
Thus, it depicts the pattern of the pulse sequence. 

The gap sequence is given by 

g. = a.-a. , , -« < i £ <« . 
y i 1 l-T 

The binary sequence is given by 

8(j) = Qf(j) -a(j-l), j £ cr . 

A smooth pulse sequence is periodic. If its rate is 
= ^ , p and q relatively prime, then the period 
is q clock times during which p pulses appear. 

Proof : (3(j+q) - 0(j ) = <r(j+q) - cKj+q-1) - cr(j ) + cr(j-l) 

= cr(j+q) - cKj) - a(j+q-l) + ct(j-I) 

= P - P = 0 . 

Lemma 5 ' makes explicit the number of pulses in each period, and Lemma 4 ' 
proves that if there were a shorter period than q , then p and q could 
not be relatively prime. 

QED 


at all others. 
Definition 6 : 

Definition 7 : 

Theorem 3: 
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These last results simplify the investigation of smoorh sequences 


in two ways. First, only rates in reduced form, that is, with p and q 
relatively prime, need be considered. Second, a single period suffices to 
describe any smooth sequence. This permits the use of the following con- 
venient notations for a smooth sequence of rate p/q : 


(l) 

Arrival times: 

( Sq t S j , . . . , * 

(ii) 

Gaps 



where q. = 

J i 

s i“ s i-i\ 


P- 1 

> indices modulo 

* 

and ^ 

%=q J 


i=0 



(iii) Counts (c(O) , ctC 1) , . . . f aCq-1) ) ; 

(iv) Binary (g(0 ) , 0 ( 1) , . . . , 3(q-l) , 


whe re 
and 


p(j ) = cr( j ) - ct(j-I) 
q-1 

I 3(J, = P 

j=0 


indices modulo q . 


Finally, we consider the difference-sequences between smooth and 
uniform sequences. 


Lemma 6 • The sequence Cs^-u^) has period q . 

Lemma b' ■ The sequence (u(j)-a(j)) has period q . 


Proofs : Since u^=^i, =u..+q for all i . By Lemma 5, fs.} satis- 

fies the same recursion, so the difference (s^-u^) has a period of q clock 
times (or p pulses). But since p and _ q are relatively prime, and 

(s i -u.) are integral only when i is a multiple of p . Therefore (s^u^) 
has no period less than q . 




With this knowfjgjdge and the obsl*'vation that all (s -lO are proper 
fractions with denominator p , and all (u(j)-a(j)) are proper fractions with 
denominator q , we can complete the refinement, begun in Lemmas 1 and 2, of 
the characterization of smooth sequences. 

Theorem 4 : A smooth sequence of rate ^ satisfies 

max(s -u.) - min (s. -u.) = . 

ill jii p 

max CuC j ) - cr(j)) - mm(u(j) - ct(j)) = . 

j , J q 

Proof : Since the difference (s^ -u.) periodically assumes p different values, 

all of which are proper fractions with denominator p , 


max( s 
i i 


-u. ) - minis -u. ) 
i ,11 


> 


2=1 

P 


This inequality combined with Definition 3 proves the theorem for the arrival 
times. A similar argument can be made for the count sequence. 

QED 


SUMMARY 


Starting from a definition of smooth sequences as synchronous pulse 
trains with bounded range of deviation from an ideal uniform rate, we have 
seen that the pattern of such a sequence is unique and periodic, and that the 
deviations from the ideal are quantized and less than unity. Dual descriptions 
have been given in terms of arrival times and cumulative counts, as well as in 
terms of their first differences. The deviations of smooth sequences from their 
uniform counterparts have been characterized. 

In closing we note that from the uniqueness of these sequences, which 
holds even under our rather broad definition of smoothness, it follows that smooth 
sequences minimize all reasonable measures of deviation from uniform spacing. 



